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SYNOPSIS 

R.L. SHANBHAG 
Ph.D. 

INDIAN INSTITUTE OP TECHNOLOGY , KANPUR 
NOVEMBER, 1980 

SOUND TRANSMISSION THROUGH SANDWICH STRUCTURES 
Tills thesis presents a theoretical analysis of the 
noise transmission characteristics of damped sandwich 
structures, with constrained damping layer treatment. The 
influence of the applied damping treatments on the struc- 
tural response of the composite structure to harmonic and 
random excitation is also studied. The following problems 
are analysed in this work. 

1. Sound transmission through an infinite sandwich plate. 

2. Noise transmission and structural response of a finite 
damped sandwich panel backed by a rectangular cavity. 

5. Noise transmission and structural response of a finite 
cylindrical sandwich shell backed by a cylindrical 
acoustical cavity. 

In addition to these problems, the problem of the 
sound transmission of a layered shell with unconstrained 
damping treatment is also analysed in this work. 



xxi 

Firstly, the sound transmission through an infinite 
sandwich plate with constrained damping layer treatment 
is investigated on classical lines. The equation of mo- 
tion of a finite three layered sandwich plate have been 
suitably reduced to represent the transverse vibration 
characteristics of the infinite sandwich plate. Closed 
form expressions for the sound transmission loss and the 
coincidence frequency are derived in terms of the sandwich 
core parameters, face layer properties and the incidence 
angle of the incident sound wave by using the dispersion 
relation between wave number and frequency. The inter- 
action between the structure and the incident, reflected 
and radiated sound fields has been considered in the ana- 
lysis. In the case of the sandwich panel, unlike the 
homogeneous panel the coincidence frequency is complex 
valued because, the structural wave motion is not a freely 
propagating wave, but a damped wave with a complex wave 
number , 

It has been shown by the analysis that, significant 
improvements in the sound insulation characteristics of 
the sandwich panel as compared to the homogeneous panel 
can be achieved, due to the shift of coincidence frequen- 
cies to higher frequency ranges and due to positive trans- 
mission loss values corresponding to the coincidence fre- 
quency. It is also revealed by the analysis that the 



xxii 


core shear parameter is the most sensitive parameter with 
respect to sound transmission and that maximum transmis- 
sion loss at coincidence can he obtained corresponding 
to an optimum value of the core shear parameter. 

Noise transmission and structural response of a 
finite sandwich panel hacked by a rectangular cavity is 
analysed next. The external pressure excitations include 
both harmonic and stationary random excitations. The 
short panel edges are assumed to be simply supported while 
the other two parallel edges are assumed to be elastical- 
ly supported. These boundary conditions represent realis- 
tic boundary conditions typical of aircraft panels sup- 
ported by orthogonal rows of stringers and frames. The 
frequency range of interest has been restricted to 1000 Hz 

An acoustoelastic formulation of the problem is made, 
where the coupling between the structural motion and the 
internal cavity pressure is fully taken into account. The 
resulting set of matrix differential equations in the 
pressure coefficients and the structural generalized co- 
ordinates is solved by a direct matrix inversion scheme 
which takes advantage of the diagonal nature of some 
of the matrices in the analysis. A 'forced damped nor- 
mal mode' analysis is made for the structural response 
in view of arbitrary damping and arbitrary boundary con- 
ditions. The natural frequencies, the associated loss 



xxiii 


factors and the damped normal modes are obtained by an 

i 

iterative interpolation method. The results are compared 
with an approximate formulation which considers only a 
partial coupling between the internal sound field and the 
structural motion. 

An exhaustive parametric study of the problem with 
respect to the sandwich core parameters is undertaken. 

It has been shown by the study, that, 

i) the boundary conditions of the panel influence the 
noise reduction inside the enclosure only in the 
low frequencies, mainly corresponding to the fun- 
damental structural resonant frequency; 

ii) the applied damping treatment improves the noise 
transmission characteristics of the sandwich plate 
significantly in the lower order structural reso- 
nant f re quenci es ; 

iii) the damping treatment has little effect on the 
noise reduction near the cavity resonant frequen- 
cies; 

iv) the structural response is reduced considerably 
due to the applied damping treatment. 

The results which are presented as noise reduction 
and structural response curves with respect to frequency, 
indicate the coupling action between the cavity and the 


structure 
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It appears from a literature survey that little or 
no work is available on the noise transmission t hr ough 
finite sandwich shells. In the thesis, the noise trans- 
mission and response characteristics of a layered shell 
with unconstrained damping treatment and a sandwich shell 
with constrained damping treatment are formulated similar 
to that of the sandwich plate taking into consideration 
the interaction of internal sound pressure field and the 
structural motion. Only the axi— symmetric modes of vibra- 
tion of the shell and simply supported end conditions are 
considered in the analysis. The problem is solved by a 
similar approach adopted for the sandwich plate. It has 
been shown by the study that the damping treatment is not 
quite effective in controlling the noise transmission and 
the structural response of the sandwich shell, because 
of low modal loss factors. 

In conclusion, it may be said, that the main contribu- 
tion of the thesis is in formulating the problem of noise 
transmission and structural response of sandwich structures 
in a general coupled acoustoelastic frame work. The the- 
sis proposes a simple matrix inversion procedure for the 
solution of the acoustoelastic problem and brings out 
clearly the use of forced damped normal modes in response 
analysis of highly damped structures . 



CHAPTER 1 


INTRODUCTION 

1.1 GENERAL INTRODUCTION , 

The interaction between sound, fields and structural 
vibration is of considerable importance in many areas of 
engineering interest, especially in aerospace engineering, 
underwater acoustics, building acoustics, and machinery 
noise control. Damage to structural parts caused by acous- 
tic fatigue, as a consequence of jet noise and boundary 
layer excitation, the transmission of noise into the inerior 
of, aircrafts causing annoyance to passengers, the vibration 
and malfunctioning of sensitive and sophisticated instrument 
packages in a vibro-acoustic environment are some of the 
problems to be considered in the design of aerospace struc- 
tural systems. In such problems, an analysis of the vibro- 
acoustic response and the sound radiation characteristics 
of structures are of great value to the designer allowing 
him to select appropriate structural parameters, geometiy, 
damping coatings and treatments to suppress structural 
vibration and control the associated radiation of sound. 

Interior noise in many transportation vehicles is 
dominated by low and intermediate frequency components. 
Especially, in propeller driven aircrafts the cabin noise 
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has maximum intensities in the low frequency range (0 - 
1000 Hz) corresponding to the blade passage frequencies. 
Acoustic absorption materials used in aircraft construc- 
tions are not quite effective in controlling interior 
noise at low frequencies [l]. Interior cabin noise at 
these frequencies, is strongly dependent on the structural 
resonances. Hence, noise attenuation in the low frequency 
ranges can best be achieved by controlling the dynamic and 
damping characteristics of structural components that radiate 
sound. 

The vibration response of structural elements, can be 
reduced significantly at resonances, by improving their 
damping characteristics. Applied damping treatments, where 
a layer of a high damping material is added to the base 
structure, have been very popular in this respect. Both 
unconstrained and constrained viscoelastic layers are used 
for the purpose. In the unconstrained damping treatment, 
the damping of the composite is achieved by dissipation 
of the vibrational energy in the extensional deformations 
of the viscoelastic layer. In the constrained damping 
layer treatment, where the viscoelastic layer of relatively 
low rigidity is sandwiched between two layers of the base 
material, the damping is mainly due to the energy dissipa- 
tion in the shear deformations of the viscoelastic core. 
Because of the reduction in vibration levels, consequent 
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reduction in the radiant acoustic energy should be possible 
with the applied damping treatments . 

Available literature on noise transmitting character- 
istics of sandwich panels with constrained viscoelastic 
damping layer is very scant. Vaicaitis [2] has considered 
the sound transmission through a viscoelastic sandwich 
panel into a rectangular cavity, with simply supported pa- 
nel edges. Recently, Reddy [3] considered the noise trans- 
mission characteristics of finite panels held in a rigid 
baffle, with both unconstrained and constrained damping 
treatments. But the emphasis of his work was on the deter- 
mination of structural response to high intensity acoustic 
excitation and was mainly experimental. The noise trans- 
mission analysis was through an approximate formulation of 
sound power calculations which were compared with his ex- 
perimental results. Moreover, in these works, the coupling 
between the acoustic field and the structural motion had 
not been treated adequately. 

Aircraft panels are usually supported by orthogonal 
rows of stringers and frames. The edge conditions obtain- 
ing in such cases cannot be treated as simple sup ports. 

A realistic formulation of the structural response and 
noise transmission problems should consider appropriate 
boundary conditions in view of the shift in the struc- 
tural resonant frequencies. In the case of panels with 

m 
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additive damping treatment, and arbitrary boundary conditions 
at the edges, determining the resonant frequencies of the 
composite structure and the associated loss factors becomes 
a difficult proposition. Moreover in such cases, the clas- 
sical normal mode analysis cannot be applied for response 
analysis. 

Statistical energy analysis (SEA), which has been suc- 
cessfully employed in structural response prediction and 
noise transmission study of homogeneous elastic structures 
in the higher frequency ranges, cannot be of much use in ana- 
lysing the low frequency noise transmission characteristics 
of viscoelastic sandwich structures. It is mainly because 
of two reasons. Firstly, in the low frequency ranges of 
interest, the modal density of the structure is small, 
making the response averages over frequency bands less 
accurate. Secondly, the assumption of weak coupling bet- 
ween modes in different frequency bands, which is widely 
adopted in SEA calculations is not quite valid for the 
sandwich panel, due to the large modal damping induced by 
the constrained damping layer, with significant intermodal 
coupling. 

1.2 SCOPE AND ORGANISATION OE THESIS 

This thesis presents a theoretical analysis of the 
noise transmission characteristics and structural response 
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of sandwich, structures with constrained damping layer treat- 
ment. The frequency range of interest in the analysis is 
restricted upto 1000 Hz in the case of finite structural sys- 
tems. However, in the study of noise transmission through 
sandwich panels of infinite extent, the restriction on 
frequency is governed by the orders of the thickness of 
the sandwich considered. The effect of the coupling bet- 
ween the sound pressure and structural motion has been con- 
sidered in. the analysis. 

The analysis admits arbitrary boundary conditions at 
the edges of the sandwich panel. These are realistic 
boundary conditions, that would obtain in sheet-stringer 
constructions commonly found in aircraft structures. Since, 
classical normal mode analyses with sinusoidal modes is not 
applicable in the case of arbitrary boundary conditions and 
arbitrary damping, a ’forced damped normal mode* analysis 
is carried out for response and transmission loss calcula- 
tions. Though the method has been outlined by Mead [4], 
perhaps, it is for the first time in this work, that a 
damped normal mode analyses is systematically made use of 
in solving a practical problem of the sound transmission 
through sandwich structures. 

The transmission loss of an infinite sandwich panel 
is considered first on classical lines. Then, the noise 
transmission through the sandwich panel into a rectangular 
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cavity is analysed. The acoustic cavity-structure inter- 
action analysis is presented in a general matrix format on 
the lines of Dowell [5]. It results in a matrix eigen 
value problem with gyroscoping coupling, characteristic of 
the dynamics of spinning structural systems. However, as 
the main interest of the present work is to determine the 
response and sound transmission characteristics of the 
sandwich panel, a direct matrix inversion scheme is adopted 
for the solution. 

There seems to be relatively very few studies available 
in the literature on sound transmission through curved struc- 
tures. Even the few works that are available are on sound 
transmission through homogeneous shells. There seems to 
be no literature on the sound transmission through sandwich 
shells. This, thesis presents a study on the vibration res- 
ponse and noise transmission of cylindrical shells with 
both constrained and unconstrained viscoelastic damping 
treatments. However, the analysis is restricted to the 
axi-symmetric vibrations of the shell and to simply sup- 
ported boundary conditions. 

In this work, two basic assumptions have been made 
regarding the properties of the layers. They are: 
i) The variation of the core loss factor of the visco- 
elastic damping material with frequency is assumed 
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to be -insignificant in the frequency range of 
interest and the viscoelastic nature of the 
material is expressed through the complex shear 
modulus G = G*(l+ip). 

ii) The base material to which the damping treatment 
is applied is elastic and nondissipative. 

The organisation of the thesis is as follows. A 
literature review of the free and forced vibration res- 
ponse of sandwich structures with constrained damping 
layer treatment is presented in the next section of this 
chapter. The review also covers the noise transmission 
studies through flexible structures, pertinent to this 
work. 


In Chapter 2, the equations of motion of a three 
layer sandwich panel is presented [6]. The method of 
determining the natural frequencies and associated modal 
loss factors for different boundary conditions is outlined. 
Expressions for forced damped normal modes, are derived 
for further analysis. 

Chapter 3 , considers the sound transmission through 
an infinite sandwich panel, for a travelling harmonic 
pressure wave excitation. The interaction between the 
sound wave and the structural motion is fully taken into 
account by considering the radiated sound on either sides 
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of the panel and the scattered sound field into the analy- 
sis. Reducing the equation of motion presented in Chapter 2 
as applicable to the infinite sandwich, the dispersion 
relation and the expression for the complex coincidence 
frequency are derived. The transmission loss (Ti) is 
obtained as a closed form expression in terms of the sand- 
wich core parameters. The effects of the core parameters 
on the TL of the infinite sandwich are presented in the 
form of graphs, T1 versus frequency. 

In Chapter 4, the structural response and noise trans- 
mission characteristics through an elastically supported 
sandwich panel into a closed rectangular cavity are analysed. 
Initially, the analysis neglects the effect of cavity pres- 
sure on structural motion and solves the problem by conver- 
ting a homogeneous differential equation with nonhomogene- 
ous boundary conditions into a nonhomogeneous differential 
equation with homogeneous boundary conditions [2], The 
solution is effected through a forced damped normal mode 
analysis. Subsequently, the coupling between the cavity 
sound pressure and the motion of the sandwich panel has 
been fully considered in an 'acoustoelastic’ formulation 
of the problem [5]. As observed earlier, a direct matrix 
inversion scheme, which takes advantage of the diagonal 
nature of some of the matrices in the acoustoelastic for- 
mulation is adopted in a partition analysis for the 
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structural response and sound transmission calculations. 

This avoids the necessity of carrying out a complicated 
eigenvalue analysis of matrices, which have complex ele- 
ments, due to the presence of the damping layer and makes 
the computational scheme much simpler. Both harmonic and 
stationary random external pressure excitations are con- 
sidered in the analysis. 

Chapter 5, deals with the noise transmission through 
finite, thin circular cylindrical shells with additive 
damping treatments. In this case, both unconstrained and 
constrained layer treatments are considered. The structural 
response including the effect of cavity is computed in the 
axisymmetric mode of vibration of the shell. Again, both 
harmonic and stationary random pressure excitations are 
considered. 

General conclusions of the present work and scope 
for further extensions of this work are presented in 
Chapter 6. 

1.3 LITERATURE REVIEW 

In this section, a review, scanning almost three 
decades of research carried out in the areas of free and 
forced vibrations of sandwich structures with additive 
damping treatments and the noise transmission characteris- 
tics of flexible structural systems is presented. 
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1.3.1 Vibration of Damped Sandwich Structures 

The free vibration characteristics of beams and plates 
with unconstrained and constrained damping treatments have 
been investigated extensively by research workers. Such 
studies include both theoretical analyses and experimental 
works and involve the derivation of the equations of motion 
of the composite structure and determination of the natural 
frequencies and the associated modal loss factors. The 
damping in these structures is obtained due to high energy 
losses mainly in the extensional mode of deformation of 
the viscoelastic layer in the case of unconstrained damping 
treatments, and in the shear deformation of the damping 
layer in the case of constrained sandwich constructions. 

Oberst [7] was the first to consider the unconstrained 
damping treatment applied to a plate. He showed that the 
damping of the composite structure depends on the stiffness, 
thickness and loss factor of the damping layer. Eerwin [8], 
presented a general analysis of damping induced by the shear 
motion of a viscoelastic layer constrained between the plate 
to be damped and a constraining layer. He showed that the 
damping factor depends on the wavelength of bending waves 
in the damped plate, the damping layer and the constraining 
layer. Ross and Eerwin [9]» developed an analysis consider- 
ing both extensional and shear damping for two layers applied 
to a plate of infinite extent. Ross, Ungar and Eerwin [10] 



11 


presented a general analysis of damping mechanism due to 
a number of constraining layers. They showed that a multi- 
layer treatment having the same weight of viscoelastic ma- 
terial has almost equal amounts of damping as in the case 
of a single constrained layer. This theory is valid only 
for simply supported conditions for finite beams. 

Using a variational approach and an energy formula- 
tion Yu [11] derived general equations of motion of sandwich 
plates and shells. Mead [6] analysed the forced flexural 
vibrations of a simply supported sandwich plate which vibra- 
ted in sinusoidal standing waves . He derived a general 
expression for the flexural stiffness. 

Heretofore listed works have assumed sinusoidal modes 
of vibration, thus restricting the theoiy to simply suppor- 
ted conditions or infinite structures. A successful attempt 
to broaden Ungar-Kerwin* s theory, for beams with finite 
length and arbitrary boundary conditions was made by 
DiTaranto [12] . He presented a general theory for free 
vibrations of sandwich beams with arbitrary boundary con- 
ditions. Complex natural frequencies w* are obtained in 
p o 

the form m * = a) (l+ip) where p is the modal loss fac- 

tor. He showed that the modal loss factor bears a rela- 

2 

tionship with w , which is independant of the end con- 

2 

ditions. However, the values of w depend on boundary 
conditions, and as such the numerical values of p need to 
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be determined separately for each boundary condition. He 
derived a sixth order differential equation in terms of 
the axial motion of the beam . 

Mead and Markus [13], showed that the modes and fre- 
quencies discussed by DiTaranto constitute a special class 
of resonant frequencies and forced modes, which exist in 
the presence of an external transverse loading, proportional 
to the local inertia loading at all points along the beam. 
Mead [14] gave a general theory of the existence of ’forced 
damped normal modes’ in structures having arbitrary damping. 
These modes can be used analogous to classical normal mode 
analysis, in steady state and stationary random response 
analysis of hysterically damped . structures. 

Yan and Dowell [15], derived a simplified equation of 
motion for the transverse displacement of sandwich beams 
and plates with constrained damping treatments. The ana- 
lytical results of natural frequencies and loss factors 
obtained by them agreed closely with that of Mead and 
Markus [16] . Mead and Markus [17,18] , in further extension 
of their earlier work determined the natural frequencies 
and modal loss factors of cantilever and encastered sand- 
wich beams. Oravsky, Markus and Simkova [19] presented 
an approximate method of calculating the modal loss fac- 
tors of damped sandwich structures by a perturbation ana- 
lysis. The theories presented by Mead and Markus [13] 
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and Yan and Dowell [15], have been found to be in close 
agreement with experimental results [20,21,22], 

The theory of vibration of sandwich beams was exten- 
ded to the case of sandwich plates by DiTaranto and McGraw 
[23], Mead [24] considered the damping properties of elas- 
tically s upported sandwich panels with a viscoelastic core, 
in terms of symmetric and antisymmetric modes. He used a 
linear 'interpolation scheme in an iterative analysis to 

determine the complex natural frequencies. Rao and Nakra 

% 

C^5] included longitudinal, transverse and rotary inertia 
effects in their analysis of flexural vibrations of 
sandwich beams and plates. Grootenhuis [26], analysed a 
multilayer damping treatment with viscoelastic materials 
and showed that unconstrained damping treatment could be 
as effective as constrained treatment compared on an equal 
weight basis. Mead [27] analysed the damping characteris- 
tics of periodic damped sandwich plates by a wave propaga- 
tion approach. Du, Killian and Everstine [28,29] used a 
finite element method in the vibration analysis <3 f three 
layered damped sandwich structures. They took into con- 
sideration the effect of frequency and temperature changes 
on the loss factor of the viscoelastic material. They 
also compared their theoretical results with experiment. 

The vibration analysis of curved members with layered 
damping treatments is not as extensive as for beams or plates. 
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Initial research on vibration problems of sandwich cylin- 
drical shells was done by Yu [30]. This work was confined 
to sandwich shells with elastic layer alone. Later, Yu[ll] 
extended his work to include layers with viscoelastic pro- 
perty by modifying the equations of motion. The shear 
damping mechanism in sandwich cylindrical shells was in- 
vestigated by Yu and Ren [31] . Jones and Salerrlj [32] also 
considered the effect of structural damping in the forced 
vibration of sandwich cylindrical shells. They showed that 
the shear damping in the core layer yields a . much lower 
loss factor for the flexural modes of sandwich shells than 
for that of the corresponding sandwich plate. ! These works 
are again applicable only to the case of simply supported 
end conditions or shells of infinite length. 

DiTaranto [33] analysed the free and forced vibrations 
of three layered sandwich rings. He used the variational 
method to derive the equations of motion and obtained ana- 
lytical expressions for the structural response consider- 
ing the circumferential bending modes. DiTaranto, Lu 
and Douglas [34] presented an approximate analysis for the 
free and forced vibrations of sandwich rings with segmen- 
ted constrained damping treatment. Lu, Douglas and Thomas 
[35] obtained the forced response of damped sandwich rings 
by a mechanical impedance calculation. They included the 
temperature and frequency effects on the core loss factor 
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in their analysis, and found good agreement with experi- 
mental results . 

Pan [36] investigated the vibrations of a damped sand- 
wich cylindrical shell in the axisymnetric mode. He follow- 
ed an analysis similar to DiTaranto [l2] and outlined a 
procedure for calculating natural frequencies and loss 
factors for finite shells with arbitrary boundary conditions. 
Markus [37] derived the equations of motion in the axisym— 
metric modes of vibration of a cylindrical shell with un- 
const rainod damping layer treatment in three configurations 
of inside, outside and two-side coatings. 

1*3.2 Sound Transmission Through Structures 

Structural response to acoustic excitation and sound 
radiation by vibrating structures and the general problem 
of sound structure interaction have been investigated by 
many author® , The different approaches to the problem can 
be found in the works of Cremer, Heckl and Ungar [38], 

Junger and Felt [39] » Lyon and Smith [40] and Lyon [41]. 

The review, here will be mainly confined to the noise 
transmission characteristics of flexible structures. 

Beranek [42], London [43], Peshback [44] and Cremer 
[45], investigated the sound transmission through flexible 
walls of infinite extent for plane harmonic wave incidence. 
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The phenomenon of coincidence transmission was presented, 
at which the wall became transparent to the incident sound 
for particular combinations of frequency and incidence an- 
gle. This is due to the matching of the trace wave length 
of the propagating sound with the free flexural wavelength 
of the plate. Subsequent works on noise transmission through 
infinite walls, divided the frequency ranges of transmission 
into three zones, namely below coincidence frequency, at or 
near coincidence and above coincidence frequency. 

The transmission of sound through double wall con- 
structions in the field of architectural acoustics was in- 
vestigated by Beranekf 46 ] , London [47], Mullohand [48], 
essentially following the lines of Beranek [42] . 

Noise transmission and structural response of cavity 
backed finite panels have received the attention of many 
research workers. Pretlove [49], using a normal mode ana- 
lysis determined the response of a cavity backed panel to 
harmonic excitation. Dowell [50], computed the response ;■>. 
and noise transmission of cavity backed finite panels to 
boundary layer turbulence. He used a numerical integra- 
tion procedure in the time domain after simulating the 
generalized forces. Dowell and Voss [51] analysed the 
response of a cavity backed finite panel clamped on all 
its edges, by a simplified modal approach, neglecting the 
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coupling between the structural modes and the structural 
acoustic coupling. Bhatt achar ya and Crocker [52], Guy 
and Bhatt achar y a [53] using a combination of Laplace trans- 
form in the time domain and Fourier transform in the Space 
domain computed the response and noise transmission of a 
finite panel backed by a rectangular cavity. They found 
good agreement between their theoretical results and experi- 
mental results . They also gave a graphical procedure to 
determine the natural frequencies of the cavity structure 
system, and explained the negative transmission loss that 
occur at frequencies close to the cavity resonances. The 
existence of coincidence like phenomenon on a cavity backed 
finite panel system was shown by Guy, Bhatt acharya and 
Crocker [54]. This happens at frequencies when the stand- 
ing wave resonant wave length of the panel matches with 
the standing wave resonant wavelength of the cavity. Unlike 
in the case of the infinite panel they showed that the coin- 
cidence frequencies of the finite panel are independent 
of the incidence angle. 

Dowell [5], and Dowell, Gorman and Smith [55] gave a 
general theory of the interaction between the internal 
sound pressure field and the elastic flexible wall of an 
enclosure and coined the term 1 acoustoelasticity’ for the 
study of such problems. They adopted a series expansion 
in terns of the structural normal modes and the cavity modes 
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and reduced "the problem of defer mining the structural 
response and the sound pressure inside the enclosure 
to a matrix eigenvalue problem with gyroscopic coupling. 
They also gave approximate results for a sinusoidal ex- 
citation, corresponding to the exciting frequency being 
close to a structural or cavity frequency, or to both. 
Recently Guy [56] gave a general analysis of the response 
of a cavity backed panel to external excitation using a 
multimodal analysis. 

The method of statistical energy analysis has also 
been used in the analysis of noise transmission through 
finite walls backed by a cavity. A balance of power flows 
into and away from the different modal systems was adopted 
in the statistical energy analysis models. Lyon [57], 
Eichler [58], White and Powell [59], Crocker and Price [60, 
61], are some of the important contributors in this respect. 
They used the formulae derived by Maidanik [62,63], Lyon 
and Smith [40] for the important SEA parameters like radia- 
tion resistances, modal densities and coupling loss fac- 
tors . 

i 

Pope [64] .and Pope and Wilby [65] investigated the 
sound transmission into closed enclosures using structu- 
ral and acoustic Green’s functions and modal analysis. 
Conditions in terms of structural and cavity configurations 
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and structural and cavity damping are established for the 
resonant and non— resonant transmission of sound into the 
enclosure. They discussed the results in relation to the 
statistical energy analysis calculations. 

McDonald, Vaicaitis and Myers [66] presented an ana- 
lytical model for the noise transmission through rectangu- 
lar plates into closed cavities using a normal mode analy- 
sis* Vaicaitis and Slazak [67] recently investigated the 
noise transmission through stiffened plates into rectangu- 
lar enclosures by a similar analysis and using transfer 
matrix methods. 

Graggs [ 68 , 69 ] considered a variational formulation 
for the structure- cavity system in the context of finite 
element representations for the coupled acousto elastic 
problem. 

The suitability of sandwich structures in improving 
the noise insulation characteristics of homogeneous panels 
was investigated by Kurtze and Watters [70]. 1 high trans- 

mission loss was sought to be achieved through a shift of 
the coincidence frequency towards higher frequency ranges 
by a choice of a thick core, whose elastic modulus is much 
less than that of the skin layers . Ford, Lord and Walker 
[71] investigated the sound transmission characteristics 
of sandwich con struct ions with a rigid polyurethane foam core. 
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They estimated the natural frequencies of the composite 
structure using an energy formulation. Their experiments 
revealed a low transmission loss due to resonance in the 
d ilat ional mode. They suggested a scheme to increase the 
transmission loss corresponding to the 'mass law* by an 
optimal choice of the core stiffness. 

Recently ,Smolenski and Erokosky [72] investigated the 
sound transmission characteristics of soft core sandwich 
panels and stressed the importance of dilat ional mode sound 
transmission. They used a Rayleigh-Ritz procedure to cal- 
culate the natural frequencies and determined the transmis- 
sion loss by an experimental procedure. Dym, Ventres and 
Lang [ 15 ], calculated the sound transmission loss of sand- 
wich panels by a theoretical analysis. They included, shear 
deformation and rotary inertia of the thick core and d na- 
tional deformation in the core axial displacement. Using 
Hamilton's equations they derived the equations of motion 
of the sandwich in symmetric and antisymmetric modes and 
calculated the transmission loss by an impedance approach. 
They also compared their results with the experimental 
results of Smolensk! and Erokosky [72]. G-uyader and 
leseueur [ 74 ] combined a statistical energy analysis and 
a normal mode approach in calculating the transmission loss 
of orthotropic multilayered plates. Madigosky and Fiorito 
[75] calculated the transmission loss as a function of 
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frequency, thickness and incident angle for single and 
multiple viscoelastic plates using a transfer matrix 
method modified to include complex wave velocities. 

Reddy [3] calculated the transmission loss of three 
layered sandwich panels, with constrained viscoelastic 
damping treatment based on a simple transmission coeffi- 
cient formula and by conversion of sound power radiated 
on the transmitting side of the panel into sound pressure 
levels . 

In all the cases of sound transmission studies through 
sandwich panels, the panels were either of infinite extent 
or confined in an infinite baffle with simple support con- 
ditions. Vaicaitis [2], considered the noise transmission 
through a sandwich plate with a viscoelastic constrained 
layer into a closed rectangular cavity. He confined his 
analysis to the frequency range of 0 - 1000 Hz and reported 
significant noise reduction. A spatially uniform and sta- 
tionary random pressure loading was considered in his ana- 
lysis . 

There seems relatively few studies on sound transmis- 
sion through curved structures. Junger and Smith [76] gave 
a semiquant itative discussion of noise transmission through 
curved structures. Smith [77], improved upon his original 
work by calculating the ratio of the absorbed power to 
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incident power of a thin cylindrical shell. Cremer and 
Heckl [78], had also considered the transmission of an 
acoustic wave through an infinite, homogeneous, isotropic 
cylinder. Their results show that the transmission proper- 
ties of the shell to an incoming single frequency wave de- 
pends on frequency and the angle between cylinder axis and 
wave normal. Furthermore, it was determined that the two 
most important characteristics of the problem were the ring 
frequency and the coincidence frequency. 

White [79] considered the sound transmission through 
finite closed cylindrical shells using a statistical energy 
analysis. He used the radiation resistance calculations of 
Maidanik and Manning [80] and the modal density formulae de- 
rived by Heckl [81] in his analysis. Koval [82,83,84] in 
a series of papers extended the work of Smith [77] on sound 
transmission through cylindrical shells by including the 
effects of external flow, reflection and scattering of in- 
cident wave, internal cabin pressurization, internal cavity 
resonances and structural damping. 

As observed earlier, there seems to be no work on 
sound transmission through sandwich shells. This thesis 
presents an investigation, on the sound transmission and 
structural response to both harmonic and random pressure 
excitation of circular cylindrical shells with unconstrained 
and constrained viscoelastic damping treatment vibrating in 
axisymmetric modes. 



CHAPTER 2 


DYNAMICS OP DAMPED SANDWICH PANEL 

The response analysis and noise transmission studies 
of the damped sandwich panel require its vibration and damp- 
ing characteristics. Different, but related analyses of the 
vibration of sandwich panels with constrained damping layer 
treatment are available in the literature [6,15,23]. Mead 
[6], derived the equations of motion of a three layered 

sandwich panel with a viscoelastic core, based on the shear 

■ r* 

damping mechanism of the core. His theory has been found 
to be in excellent agreement with experimental results [20, 
21,22] in the usual ranges of frequency interest and widely 
adopted. In the present work the- equations of motion deri- 
ved by Mead are adopted in the study of structural response 
and noise transmission of sandwich panels. They are presen- 
ted in this chapter for the sake of continuity and under- 
standing of the subsequent chapters. The steps involved 
in the determination of the natural frequencies and the 
associated modal loss factors and the response analysis 
using damped normal modal analysis are also briefly outlined. 

2.1 EQUATIONS OP MOTION 

The equations of motion of a uniform three layer sand- 
wich panel under transverse pressure excitation was derived 
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by Mead [6], under the following assumptions. 

i) Direct stresses in the soft core are much smaller 
than the direct stresses in the face plates and 
hence neglected. 

ii) There is no significant direct strain in the core 
perpendicular to the plane of the plate. Therefore, 
the plate motion in the transverse direction w is 
the same for all the three layers. 

iii) There is no significant shear strain in either face- 
plate, in planes perpendicular to the plane of the 
plate* 

The equations of motion, neglecting the inertia forces 
in the x and y direction are, (Refer to Figure 2*1) 
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where , 

B = Young* s modulus of the plate 
v = Poisson’s ratio of the plate 

D t = total flexural rigidity of both the face plates 
= [2Eh^/l2 (l-v 2 )j 

& = G- * (l+ip ) , complex shear modulus of the core 

P = core loss factor 

M- = mass per unit area of the whole sandwich panel 
»"k ) = local resultant transverse loading on the 
sandwich plate 

V 2 „ _aj + aL v 4 a * 3 4 E 4 

dy 2x 4 Cx 2 6y 2 2y 4 

u = midplane displacement in the x-direction of the 
bottom face plate 

— midplane displacement in the y— direction of the 
bottom face plate. 

In the above equations it has been assumed that the 
top and bottom face plates are of the same material pnd 
thickness . 

It is convenient for subsequent analysis, to express 
the equations of motion in terms of the transverse dis- 
placement w, by eliminating u and v from equations (2.1) 
to (2.3). This is done in following mann er. 
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Equations (2.1) to (2.3) can be rewritten as, 
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where , 
Y 

S' = 


3 (l+l^/t^ ) 2 
2 G- (1- v 2 ) 


is the geometric parameter 
= g(l+ip) is the complex shear parameter. 


Differentiating equation (2.4) and (2.5) with respect 
to x and y respectively, and adding the resulting equations 
one obtains , 
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Differentiating equation (2.6) twice with respect to 
x and y separately, and adding the resulting equations, 
the following equation is obtained. 



28 


- S' Y V 4w + 


(V 2 w) + — ( ^- h2 - (^f 

B t St h£ Sx 5 


+ -£% + 
3xdy 


S^V , 8^-y \ 

2T h ; = 

SxSy Sy^ 


2 

VP 


(2.8) 


where, 


V^ is the differential operator 
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Using equations (2.6) and (2.7) to eliminate the u and 
v terms from equation (2.8) one gets. 
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It has heen shown hy Mead [13,14], that separable 
harmonic motions for w with characteristic frequencies 

are possible, when the structure is excited by a particu- 

a ■ 

lar form of forcing functions called the ’damped normsl 
loading* . The corresponding modes in which the plate 
vibrates are called the forced damped normal modes. The 
damped normal loading for the sandwich panel is harmonic 
and is in phase with the local velocity but proportional 
in amplitude to the local inertia force. Hence the damped 


normal transverse loadings are of the form 
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A 

p(x,y,t) = -it^p. -Mf (2.10) 

dt^ 

where, 

t) = the loss factor associated with a damped normal mode. 


Under the damped normal loading, the plate motion can 
he expressed in the form, 

w = W(x,y) e ltot (2.11) 

where, 

¥(x,y) represents the forced damped normal mode which is 
generally a complex function of x and y. 


Substituting equations (2.10) and (2.11) into equa- 
tion (2.9), one obtains 


V 6 W - g'(l+Y) V 4 w- “ 2 (l+it)) — [ ^ 2 W - g'W] = 0 (2.12) 

D t 


Equation (2.12) can be rewritten in non-dimensional form 
using non-dimensional coordinates % = x/a and ^ = y/a, 
where a is the length of the plate in the x-direetion. 

The result is 

V 6 ¥ — g( 1+Y ) V 4 W— d V 2 ¥+g dll =0 (2.13) 


In the above equation the differentiations are with 
respect to the nondimensional coordinates, - and ^ in 
the operator . 
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2 

g = g'a = g* (l+ip) is the nondimensional shear parameter 
and 

p yj p 

n 2 = - - - Xl+fo) ) p. a _ (l+irj) is the nondimensional 

"fc 

frequency parameter. 

2.2 NATURAL FREQUENCIES , LOSS FACTORS AND DAMPED FORMAL MODES 

In the analysis, panel edges are assumed to be simply 
supported along the x-wise edges at y = 0,b. The y-wise 
edges at x = + a/2 are assumed to be elastically supported. 
These assumptions are justified on the following grounds . 

The noise transmission studied in the present work, is in 
the context of interior cabin noise of aircrafts. Typical 
aircraft panels are supported by orthogonal rows of strin- 
gers and frames. It has been found that the boundary con- 
ditions at the x-wise edges have little effect on the na- 
tural frequencies of the plate, provided the width of the 
panel between the frames (in the y-direction) is more than 
two or three times the stiff ener.. spacing(in the x-direction). 
Moreover, previous works on the vibration analysis of sheet 
stiffener combinations have consistently made these assump- 
tions [85,86,87]. The analysis can as well include the 
admission of arbitrary boundary conditions along both x and 
y—wise edges, but in such a case the computational efforts 
in determining the complex eigenvalues and the damped nor- 
mal modes would become enormous . 
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Under the above assumptions, the forced damped normal 
modes can be expressed as 

6 w 

, nita? 
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r=l b 

Substituting equation (2.14) in equation (2.13) a 

characteristic polynomial equation in A is obtained. 
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Writing u and v analogous to equation (2.11) as. 


u = U(x,y) exp(icot) and v = V(x,y )exp(i w t ) 

the damped normal modes corresponding to displacements 

u and v can be expressed as, 

' 6 

= S B rmn * > 8± n P r, * (2.16) 


V (£,¥)= 2 C exp( X 5 )eos P ¥ 

mn v ’ ' u rmn - v rmn n 


(2.17) 


Substituting equations (2.14), (2.16) and (2.17) into 
the equations of motion (2.1) and (2.2), it can be shown 


that the coefficients B___ and C 

rmn rmn 

the following form. 


related to A rmn in 


are 
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B 


rmn 


A nm 6 (h l + V . 

2a( A 2 - p^ - g ) ■ rmn 
v x rmn n ^ J 


n _ 

xian ~ 


P n s + ^ 2 ) 
2sl ( ^rmn P n ~ 


■^rmn 


(2.18) 

(2.19) 


The problem of determining the natural frequencies, 

modal loss factors and the damped normal modes reduces to 

finding the set of values A , P and that satisfy 

rmn n 

both equation (2.15) and the boundary conditions at x = +a/2. 


2.3 BOUKDABJ CONDITIONS ALONG Y-¥ISE EDGES AT x = + a/2 

The following boundary conditions along the edges at 
x = + a/2 are considered in the present analysis. 

i) Edges simply supported 

ii ) Edges clamped 

iii) Edges with elastic rotational restraint and no trans- 
verse displacement, with and without rivets 

iv) Edges with no rotation but with transverse elastic 
restraint. 

The nature of the boundary conditions in the above 
cases are given in reference [24]. 

2.3.1 Simply Supported Edges 

In this case the natural frequencies , t he modal 
loss factors and the damped normal modes can be obtained 
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in closed form. The boundary conditions are. 


w 1 SB 0 

*x = + a/2 

(zero transverse displacement) 

+ V » J 5 ) +a/2 ’ ° 

(zero axial strain in entire 
face plate) 

M x s* 2M^ + 

(zero bending moment per 


unit length along the 

y-wise edges of the plate) 

+V ^ L ±a/2 = ° 



Applying the boundary conditions (2.20 ) it can be shown 
that the damped normal modes in the case of simply supported 
conditions at x = + a/ 2 are given by 


\ n (- » = A mn sin nra( Z+ 1) s in (2.21) 


Closed form expression for the nondimensional natural fre- 
quency parameter can be obtained by substitution of equation 
(2.21) in equation (2.13). 


The result is 


= SS < 1+1 W “ 


(m 2 n 2 +P 2 )+g(l+Y )(m 2 rc 2 +P 2 ) 2 
m 2 % 2 + P* + g 


( 2 . 22 ) 
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The real part of the above equation gives the natural 
frequency parameter and the ratio of the imaginary part to 
the real part gives the modal loss factor, 

2.3.2 damned Edges 

The boundary conditions in this case are: 

(zero transverse 
displacement ) 

(zero slope) 

(zero shear strain in 
the core in x~z plane) 

Applying the above boundary conditions and using 
equations (2.14) and (2.16), the natural frequencies, loss 
factors, and the damped normal modes can be obtained, by 
solving the resulting determinant al equation in conjunction 
with equation (2.15). The natural frequencies are those 
which satisfy the equation (2.15) and make the determinant 
of the coefficient matrix of A rmn zero. 

2.3.3 Edges with Elastic Rotational Restraint and No 
Transverse Displacement 

Such boundary conditions for the panel will arise when 
the edge stiff eners have finite rotational stiffness and 
very large transverse stiffness. 


w | = 0 

x= +a/2 


3w 

ax 


I 


= 0 


x= ±a/2 


**2 




= 0 


x = ±a/2 
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a) For the ease of unriveted edges the boundary con- 
ditions become 


w I =0 ( zero transverse displacement) 

x = ± a/2 

It will be assumed that at an unriveted edge there 
are no mid -plane direct stress cr x in either face plate. 
This implies that 



x = +a/2 


The last two boundary conditions are due to the mo- 
ment equilibrium between the stiffener and the plate and 
are expressed as, 


v 


A 

~ 

Sxdy c 


E P 

3 


a 5 w 

A 

2x3y 


I 

x = +a/2 



+a/2 


where G J is the St. Tenant torsional stiffness of the 
s 

stiffener. Pis the Wagner-Kappus torsion bending con- 
stant, E 0 is the stiffener Young’s modulus, and M_ is the 
moment per unit length along the y-wise edge of the plate 
to which the stiffener is attached. 

b) For the case of riveted edges, the shear strain in the 
core is completely prevented. In this case the six boun- 
dary conditions are: 
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w | = 0 

x * + a/2 


(zero transverse displacement) 


|-[-2u 4 -(h 1 +h 2 )||] | 

- x = ± a/2 


=0 


(zero shear strain in the 
core in the x-z plane) 


and 


G J 
s 




dxay 


2 - E s r 


.IJL 


Qx5y 


x= +a/2 


(moment equilibrium bet- 
ween plate and stiffener) 



±a/2 


As for the clamped edges, the natural frequencies, 
loss factors and damped normal modes are obtained after 
solution of a determinantal equation in conjunction with 
the polynomial equation in As. 


2.3.4 Edges with Ho Rotation but Transverse Elastic Restraint 

Such boundary conditions for the panel will arise, 
when the edge stiffeners have finite transverse stiffness 
and have large rotational stiffness. 

Zero rotation at an edge is expressed as 

—Ip | = 0 (zero rotation) 

x x = + a/2 

Two more boundary conditions are due to equilibrium 
of the shear forces between the plate and the stiffener 
and are expressed as 
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SM^ 

dx 


x 2 


eM 


SL 


3y 


+a/2 


?»4 

= (SI) 

oy 4 


x = + a/2 


where, (SI) is the flexural rigidity of the stiffener 


= -[D t (l-V ) xfw + Jh 


8xdy 2(1+ v ) 




where, M is the twisting moment per unit length along 
the y— wise edge of the plate. 


The last two boundary conditions corresponding to the 
riveted edge conditions are 


fe; t~ 2n + + i^) ll ] \ x = + 


a/2 


= 0, (zero shear strain 
in the core in 


2.4 SYMMETRIC MODES 

If the boundary conditions are identical at x = + a/2 
and if the loading is symmetric, the contributions to the 
total response will be mainly due to the symmetric modes and 
it is expedient in such cases to consider only the symmetric 
modes of vibration, which can be expressed in the form 

’W ! '* ) = S V °° 8 *rmn« 
r=l 

5 * * 

U ( f , n) = £ B sin A S 

mn v nun ran 

r=l 


sin? n f (2.23) 

sin 


(2.24) 
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W 5 . *> 


(2.25) 


In this case the coefficients B , C and A 

mm nan rmn 

are related as. 


B. 


rmn 


nan 


. S A mn * b 2 ) 

2 a (g + A 2 + P 2 ) 
ran n 

6 g n (y» y 
2 a (« + A Ln + S n> 


^rmi 


^rn 


n 


(2.26) 


(2.27) 


and the polynomial equation in ' s becomes 

{ A 2 + P 2 )^ + a: ( i+y ) ( A 2 + P 2 )-n 2 ( A 2 +P 2 + g) = 0 

v rmn r n y s ' lX ' rmn n ; ' rmn n 6/ 

(2.28) 


The effort of computing the natural frequencies is 
reduced considerably as the boundary conditions are to be 
applied only at one edge namely at x = + a/2 and the fre- 
quency determinant reduces to the order 3x3. 

The elements of the coefficient matrix for the dif- 
ferent boundary conditions and the forced damped normal 
modes are given in Appendix A. 


2.5 DETERMINATION OF RESONANT PREgJENCIES AND LOSS FACTORS. 

The solution of the determinantal equation is compli- 
cated because of the cumbersome complex arithmatic involved 
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and the inherent numerical problems associated with it. 

This difficulty can be overcome by splitting the deter- 
minant into its real and imaginary parts and setting these 
separately equal to zero. Considering the real part and 

the imaginary part of the determinant as two real functions 

2 2 

in the two real variables and y] q * , a two dimen- 

sional Newton-Raphson iteration scheme as reported by 
Narayanan and Mailik [88] can be adopted. But a direct 
iterative linear interpolation scheme as given by Mead [24] 
is used because of faster convergence. The steps in the 
procedure are as follows: 

i) A value of one of the A.J. s i s guessed, say A_^. 

For a first approximation, the imaginary part may 
be neglected. 

.,2 

ii) The first guessed value of A^ was used to find Cl 
from equation (2.15) or (2.28). Either of the cor- 
responding two equations (2.15) or (2.28) is used 

to find the approximate values of the otter two A r s - 

iii) With the three approximate values of A rmn , the 
determinant of the coefficients for the approxi- 
mate boundary conditions is evaluated. 

iv) The value of A is slightly perturbed ate steps(ii) 
and (iii) are repeated to find another value of the 
determinant. Linear interpolation is adopted to 
determine subsequent values of A which brings 
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the determinant closer to zero. When the X s 

2 

converge sufficiently, the corresponding value of Cl 

gives the frequency parameter and loss factor. 

If the initial guess of X is properly made then the 
convergence is obtained within four to five iterations. 

In sane cases the number of iterations for convergence 
increased for higher order modes. In any case the number 
of iterations for convergence did not exceed twenty. 

After getting the correct values of A s, and sub- 
stituting into equations (2.14) or equation (2.23) the dam- 
ped normal modes are obtained. The form of the damped 
normal modes in the symmetric modes for the different 
boundary conditions, are given in Appendix A. The values 
of the natural frequencies and loss factors are given in 
Chapter 4, for the typical plate dimensions and layer pro- 
perties adopted. \ 

2.6 FORCSD DAMPED NORMAL MODE RESPONSE ANALYSIS 
2.6.1 Harmonic Loading , 

Considering a harmonic external pressure 
p(x,y,0,t ) ss exp( i oj t), the damped normal modes can be 
made use of to obtain the response of the sandwich plate. 

The forced response of the plate can be expanded in 
terms of the damped normal modes as 
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w(x r y,t) 


©O CO 




y«> v (t) 


(2.29) 


Substituting equation (2.29) into equation (2.9) and 
using orthogonality properties of the damped normal modes 
[14] > it can be shown that 




yz 


p ( t ) 
%n s J 


u [ u (1+ i n ) 

nnn L in 'in 


— 2 — 

« 3 


(2.30) 


where , 

= 


. m 

b a/2 

/ ' / p(x,y,0,t)W _ (x,y) dx dy 

o —a/2 mn 


(2.31) 


and 

• b a/2 „ A 

Pan * / / /0 P < x »y> ix (2.32) 

o —a/2 . ^|J 

The p mn ' s and s can be considered as the generalized 
force and generalized mass analogous to classical normal 
mode analysis. The expressions for p mn and p ^ correspond- 
ing to the different boundary conditions are given in 
Appendix B. 


2.6.2 Random Loading 

Consider a random pressure excitation stationary 
with respect to time and homogeneous with respect to the 
spatial indexing parameters, let the cross power spectral 
density of the pressures at two locations x^,y^ 
on the panel be expressed as 
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S pp (x L' y l ; x 2» y 2» a >. } = S p 0 p 0 (a) )r xy (x 2“ x l iy 2" y l’ w } 

(2.33) 

where , (x^-x^; y 2 -y-^; u ) is in the form of a correla- 

tion function and being a function of <0 may be regarded as 
the spectrum of the correlation coefficient. 


It can be shown that [89] the cross power spectral 
density of the displacement w is given by 


S ww^ s l* y ! i *2 




111=0 




H ( 10 ) H* ( <0 ) I (( 0 ) (2.34) 

mn' ’ rs v y mnrs v J j 


where, 

*W“> 


u mn^ +it W^ - (0 ] 


is the complex 


th 


frequency response function in the mn damped normal mode, 
the asterisk denotes the complex conjugate and 
b a/2 b a/2 

We < “ >- sj /z f_f /z y *i .?! )*»< *2 .y 2 ) 


r xy^ x L ,y l i x 2» y 2 ; 05 ^ dx 1 dy 1 dx 2 dy 2 (2.35) 

Ihe auto power spectral density of the displacement 
corresponding to a point (x,y) in the panel is obtained from 
equation (2.35) by setting x^x^x and y 1= y 2 =y. Expressions 
for I mnrs (<0 ) are also derived in Appendix B for different 
boundary conditions of the plate and for a uniform random 
pressure. ( r xy ^ x i»yi» M ^ = 1-0). 


CHAPTMt 3 


SOUND TRANSMISSION BY DAMPED SANDWICH 
PANELS OF INFINITE EXTENT 

3.1 INTRODUCTION 

The interaction of an infinite homogeneous panel 
and and an airborne sound wave incident on one side of 
the panel is a classical problem in acoustics. The pheno- 
menon of coincidence, when the bending wavelength of the 
panel matches with the projected wavelength of the sound 
wave on to the panel, resulting in a high degree of coup- 
ling between the sound and structure is explained in the 
context of this problem. For grazing incidence, the fre- 
quency at which perfect matching of the free flexural wave- 
length and the sound wave length occurs is called the cri- 
tical frequency and at this frequency for grazing incidence, 
the plate radiates sound of equal intensity as the inci- 
dent sound in the transmission side of the panel. For other 
incident angles , the wave coincidence occurs at frequencies 
greater than the critical frequency depending on the inci- 
dent angle. 

In the present chapter, the sound transmission charac- 
teristics of a damped sandwich panel of infinite extent is 
investigated on classical lines. The sandwich panel consists 
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of two face plates of the same material and thickness, 
with a constrained viscoelastic damping layer in between. 

The expression for sound transmission loss in terms of the 
core parameters is derived, taking into account the inci- 
dent, reflected and radiated waves from the panel. The con- 
dition for coincidence transmission is established in terms 
of core parameters and the incident angle of the exciting 
pressure, using the dispersion relation. 

Parametric studies on transmission loss are carried 
out and the results are presented in the form of transmis- 
sion loss curves. It is shown by the analysis that the 
core shear parameter has significant effect on the noise 
transmission characteristics of the sandwich panel, and 
the location of the coincidence frequencies. The core loss 
factor and the geometric parameter also influence the trans- 
mission characteristics of the sandwich panel. Results are 
compared with an equivalent homogeneous plate of face plate 
material having the same surface mass density as the sand- 
wich. It is seen from the results that the sandwich panel 
has better sound transmission characteristics based on two 
c o tints . Firstly, the coincidence frequencies are shifted 
to higher frequency ranges for the case of the sandwich 
panel. Secondly, the coincidence transmission loss is con- 
siderably improved, due to the presence of the damping 
layer. 
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3.2 S QUITO TRANSMISSION LOSS 

Consider the sandwich panel of infinite extent, se- 
parating two semi-infinite air spaces as shown in Figure 
3.1. Let a plane harmonic pressure wave he incident on 
the panel at an angle 8 . 8 is the angle between the 

normal to the plate and the normal to the wave fronts. 

The incident pressure can be expressed as 
p^x^t) = P ± exp[ik(x sin 0 -zcos e )]e i (3.1) 

where, k is the wave number of the incident wave and 
P^ is the amplitude of the incident pressure . In the ana- 
lysis to follow, the harmonic term is consistently 

dropped from all equations as a common term. 

The resultant pressure on the incident side of the 
panel consists of the incident pressure p^, the reflected 
pressure p and the radiated pressure denoted by p . The 
reflected pressure p g corresponds to the scattered pressure 
from the panel when it is held rigid. The radiated pres- 
sure on the transmission side denoted by p^ is the trans- 
mitted pressure. 

The reflected pressure p g can be shown to be, [39] 

p s (x,z) = exp [ik (x sine + z cos e )] (3 2 . ) 

The difference between the expression for the reflected 
pressure and the incident pressure is only in the sign of 
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the z coordinate which is reversed, implying thereby that 
the wave is specularly reflected. 

As the external pressure excitation on the sandwich 
panel is independent of the y coordinate, the dynamic res- 
ponse of the plate must also he independent of y. This 
implies that the dynamic response consists of structure- 
borne waves propagating in the x direction, the wave fronts 
being parallel to the y axis. 

The plate acceleration can be written as, 
w(x) = ¥ exp (ilex sin 9 ) (3.5) 

The transmitted pressure p^ can be determined by solving 
the equation of motion to the resultant surface pressure 
loading on the panel. Since the radiated pressures on 
either side of the panel are dependent on the structural 
motion, which in turn depends on the resultant pressure 
loading, the problem of determining p^ represents a coupled 
structural -acoustic problem. 

The radiated pressure p r satisfies the two dimen- 
sional Helmholtz equation which for the acceleration dis- 
tribution represented by equation (3.3) takes the form 
a 2 « 2 ~ 

(-^s- + —5 +■ k ) p(x,z) = 0 (3.4) 

3x^ 6z 

On the panel surface, the pressure satisfies the boundary 


condition 
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3P 

= ~pW exp (ikx sin 9) j 

32 z = 0 

where p is the density of fluid medium (air). For 
waves propagating in the positive z direction, the pres- 
sure field satisfying the two dimensional Helmholtz equa- 
tion and the boundary- condition at the surface is of the 
form, 

• ipW 

p ( x, z) = exp [ik(xsin 8+ z cos 8)] (3.6) 

kcose 

The above expression represents the radiated pressure 
on the incident side. 

The transmitted pressure p^. associated with the elas- 
tic response of the plate is the pressure radiated by the 
plate into the half space z<0. Since the two panel sur- 
faces on the incident and transmitted sides have the same 
transverse motion, and the acoustic medium on either side 
of the plate being the sane 

|Ptj = |Prj (3 ' 7) 

On the plate surface , on the transmitted side, the transmi- 
tted pressure p^ takes the negative value of the radiated 
pressure on the plate surface on the incident side, since 
the plate acceleration W(x) producing a positive pressure 
on one plate surface generates a negative pressure on the 
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other surface. As the transnitt ed wave propagates in the 
negative z direction, it can he expressed as 


P t (x,z) = | p t | exp [ik(xsine- zcos 0)] 

, , , % (3.7) 

= -p r (x,z) exp (~i2kzcos 0 ) 

The resultant plate loading is now obtained as the sun of 
the incident, reflected and the radiated pressures on either 
side of the panel at the panel surface. The pressure load- 
ing on the plate surface is then, 

• % ’ ' 

p(x,0) = ~P ± “ P s - P r + P t | = -2 Pj , - 2p r (3.8) 

z = 0 


The equation of notion for the sandwich panel considered 
in the analysis, is given by equation (2.9). For the y 
independent configuration it reduces to 


5^w , / 1V \ 5 4 w ^ r 

g ( 1+X ) J + *r— [ 

.6 dx 4 r Sx^St 


8 4 w 

~T77Z 


3x 


g TJ 

8 1 ^ 


1_ 

D, 


( Hr - S’p) 


't w X~ 

where the various terns representing the sandwich core 
parameters have been explained in the previous chapter. 
From equation (3.3), for harmonic motion, the displace- 
ment can be expressed as. 


(3.9) 


yr 

w(x) = - — p exp (ikx sin 0) 


(3.10) 



50 


Substitution of equation (3.10) in equation (3.9), using 
equation (3.1) yields. 


£ r k 6 sin 6 e + g» (1+Y )Ain 4 e _ fc 2 

M 2 u 2 “ jv 


sin 2 9 


_ 


kcos 0 D, 


(k 2 sin 2 9 + g» )] 


2 p 


B. 


— (k 2 sin 2 0 + g’ ) 


(3.11) 


Solving for ¥, we get. 


¥ = 


2 P^^ cose 


^#- e a + f^ — )-x| oos e i 

** l ^ ^ r sin 9 +g' i 


(3.12) 


Hence the radiated pressures on either side of the sand- 
wich panel can be obtained frcm equations ( 3 . 6 ) and ( 3 . 7 ). 


p„(x,z) = 


2 p. esqj [ik (x sin9 + z cos 0 ) ] 


ikp 

~T~ 


D t k 4 sin 4 _e ( )-ll cose -2 


( 1 ) 


k sin e+g’ 


(3.13) 


and 

P t (x,z) 


2 P^ exp [ik (x sin 9 - z cos 0 )] 


P ( ^ 0 )^ k i 


£1 H i 

P (. » a 2 ^kWw' ) 


cos 0-2 


(3.14) 


From equations (3-14) and (3.1) the complex ratio of the 
incident and transmitted sound pressures for the sandwich 
panel can be writt en as , 
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P i 1_ l^cosSj- ^t kjsinla (l+ g 

w k sirtg-irg 1 


jjo. 


P t 2P - li 

The above equation can be written in the form 


O-l] (3.15) 


li 

H 


i + 


Zrp cose 
2 pc 


(3.16) 


analogous to the case of a homogeneous panel, where can 
be considered as the specific transmission impedance of the 
damped sandwich panel. Unlike in the case of the homogene- 
ous elastic panel, where the specific transmission impe- 
dance is always reactive, in the case of the sandwich panel 
with viscoelastic layer, the specific transmission impedance 
has both reactive and resistive parts. can be represen- 

ted as 


— Efji + 1 Xfp 


(3.17) 


where is the specific transmission resistance and 2^, is 
the specific transmission reactance. The expressions for 
and 2^ can be expressed in terms of the sandwich core 
parameters and the properties of the elastic face plates 
and the incident angle 


k 6 sin 6 e D t g p I 

= , - 2 -- T -~ — 

w(k sin e+ g ) 4-g . p 

and 


(3.18) 
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D t k 4 sin 4 e £?(k 2 sin 2 e + g + ^ 2 ) 

*[l+ -7-5 5 9 — 9 9 1 




(k 2 sin 2 e +g) 2 +g 2 p 2 


(3.19) 


where g is the real part of g* 


Defining the sound transmission coefficient t of the sand- 
wich panel as the ratio of the transmitted sound power to 
the incident sound power, one can express 


_ 1 _ 

T 




(l+ cos 9 Y" -t- ( 

P ' 




2 ^ c 2p c 


cos 0 )‘ 


(3.20) 


and the sound transmission loss of the sandwich panel can 
he expressed as 


Ti = 10 log 


( 3 . 21 ) 


It can he verified that equation (3.21) reduces to the trans- 
mission loss formula of singLe homogeneous panel, whose bend- 
ing rigidity is twice the bending rigidity of either of the 

face plates as g* tends to zero, where the specific trans- 

a 

mission impedance becomes wholly reactive. For an elas- 
tic core with {3=0 also, the specific transmission impe— 

Cl,.. 

dfnce becomes reactive at all frequencies,. 


3.3 COINCIDENCE FREQUENCY AND COINCIDENCE NOISE TRANSMISSION 

Unlike the case of the homogeneous elastic panel, for 
the damped sandwich panel, the coincidence frequency becomes 
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a complex number, as it depends on the core parameters. 
Because of the presence of the damping layer, propagation 
of free flexural waves in the sandwich panel is not pos- 
sible. However, a damped wave motion decaying with res- 
pect to the spatial coordinate is possible. It can be 
thought of as a forced damped wave corresponding to the 
forced damped normal modes of finite sandwich panels dis- 
cussed in Chapter 2. Hence, the structural wave number or 
the propagation constant becomes a complex number. Since 
coincidence is the condition of matching of the projected 
acoustic waves on to the panel and the structural wave 
forms, the coincidence frequency also becomes complex. 

The real part of the coincidence frequency represents the 
frequency of minimum transmission loss and the imaginary 
part indicates nonzero transmission loss at coincidence 
for the sandwich panel. The coincidence frequency can be 
derived for the sandwich panel by considering the disper- 
sion relation between frequency and wave number. 

Assuming a wave solution of the form 

w = W exp (ikx) (3.22) 

and substituting into the equation of motion (3.9) with 
the external loading terms zero, we get, 

-k 6 - g* (1+Y) k 4 + w 2 k 2 + ~ 2 g» =0 

*fc t 


(3.23) 
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In the foregoing equation, it should be noted that the 

wave number k is a complex quantity. Writing k =w /eg, 

where Cg is the flexural wave speed in the panel, equation 

o 

(3.2» represents a quadratic in , the solution of which 
can be written as 

2 C B ^ 2 

(a = £[ - Cg g* (l+Y)] 

T\ -O' 


±4- V([ - 4 g’( 1+Y )3 2 + ' — ) 

D+. * D. 


(3.24) 


Since the coincidence occurs when the trace acoustic 
velocity on the panel matches with the flexural wave speed, 
the coincidence condition is established by setting Cg=c/sin 0 
in equation (3.24) where c is the local acoustic velocity 
in the medium. Admitting only the positive real part for 
the complex roo"^ in equation (3.24), the coincidence fre- 
quency for different angles of incidence can be expressed 
as 

4 2 
o c ji c 

» c = * [ T~ TT «' < 1+Y) 3 

D^sin e sin q 

4 2 . , 6 
cV c p 4g'iiC 

+ i V([ J z- g'(l+Y)] + c~) 

D^sin 0 sin 0 sin @ 

(3.25) 

The real part of the above equation gives the coincidence 
frequency for the sandwich panel and the ratio of the 
Imaginary part to the real part represents a measure of 
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the rate of decay of the damped wave in the structure, 
the modulus of which nay be likened to the loss factor. 

The same expression for the coincidence frequency 
equation ( 3 . 25 ) can also be obtained by requiring the 
amplitude of the transmitted pressure p^ to be equal to 

P i* 


The minimum value of coincidence frequency is obtained 
by letting 9= tr/2 in equation (3.25), which is the critical 
frequency of the sandwich panel being the coincidence fre- 
quency for grazing incidence. 


The phase velocity of bending waves of the sandwich 
plate is obtained from equation (3.23) as 


<4 - [ k2 + i 

V (g*/k 2 +1) 


( 3 . 26 ) 


which is the dispersion relation for the sandwich panel. 
The above relation reduces to the dispersion relation of 
a homogeneous panel having flexural rigidity equal to twice 
the flexural rigidity of either of the face plates as the 
shear parameter g* tends to zero, that is, with a very soft 
core. As the thickness of the core hu, becomes very small, 
(hg - * 0) the flexural wave speed corresponds to the flexu- 
ral wave speed of a homogeneous panel having twice the 
thickness of the face plate. For practical values of the 
shear parameter and geometric parameter and relatively 
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higher frequencies, equation (3.26) can be approximated by 
4 * ~ (3.27) 

p 

which is very similar in form as the flexural wave speed of 
homogeneous panels. 


3.4 RESULTS AND DISCUSSION 

The transmission loss for the infinite sandwich panel 
has been obtained from equation (3.21), for different values 
of the core parameters, viz, the core loss factor, the shear 
parameter and the geometric parameter. The face plate con- 
sidered in the analysis is aluminium having Young 1 s modulus 
E = 7.24 x 10^ E/m^ and density p = 2770 kgm/m^ . The 
viscoelastic core considered corresponds to typical PVC 
materials, the variations in whose properties have been 
taken care of by the parametric study of the core shear 
parameter g and the loss factor p. The storage shear 
modulus of the core and its loss factor have been assumed 
to be frequency independent. The face plate is assumed to 
be elastic and nondissipative . 


The local acoustic velocity in air is taken as 
c = 330 m/sec. and the density of air is taken as 
p = 1.225 kgm/m^. 
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The thickness of the face plate considered for para- 
metric study is h^ = 1 mm. Variations in the sound trans- 
mission loss with changes in thickness of the face plate 
have also been studied. For a thin panel, such as consi- 
dered in the analysis ,equat ion (3.9) quite accurately pre- 
dicts the behaviour of the sandwich panel as the flexural 
wavelengths are many times larger than the thickness of 
the face plates (about 30 to 40 times) in the frequency 
range of interest considered [16,21]. It obviates the 
necessity of including the effects of shear deflection 
and rotary inertia. Even for flexural wavelengths which 
are about twenty times the thickness of the face plate, 
the approximations implied in equation (3.9) by neglect 
of the shear deflection and rotary inertia effects will 
not significantly impair the accuracy of the results. The 
range of frequencies for which the equation of motion will 
be valid without much error can be determined from equa- 
tion (3.26) for the face plate thickness adopted. For 
example, h^ = 1 mm the range of frequency, for a core of 
the same thickness as the face plate and g = 100 turns 
out to be 10440 Hz, for bending wavelengths thirty times 
the thickness of the face plate. The valid range of fre- 
quency will increase with increasing values of the core 
shear parameter g as is evident from equation (3.26). 

The density of the core material p is taken to 

c 

be one-fifths of the density of the face plate material 
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2 

( P c = 554 kgm/m ). in all calculations. 

The coincidence frequencies and the loss factors 
calculated from equation (5.25) are shown in Table 5.1 
for various incidence angles and the core shear parameters 
g . Prom the table it is observed that the coincidence 
frequency decreases with increasing values of the core 
shear parameter g for all angles of incidence. It is also 
seen that the coincidence frequencies decrease with increa- 
sing angles of incidence as it should be from equation( 5 .25 )* 
having minimum values for grazing incidence. In the prac- 
tical ranges of frequency coincidence conditions exist only 
for incidence angles greater than about 45 °. 

This observation is also true for homogeneous elastic 
panels where the upper limit for the coincidence frequency 
is given by the condition for Rayleigh waves [90], which 
places a lower limit on the incidence angle 9 for which 
coincidence can occur. It is also seen from the table 
that the loss factor initially increases with increase in 
the core shear parameter g , upto a certain value of g 
and then decreases for further increase in g . This is 
especially seen for incidence angles greater than 45°. 

Thus the value of g at which the maximum loss factor is 
obtainld also depends on the incidence angle q . The 
existence of an optimum value of the core shear parameter g 
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for maximum damping effectiveness as revealed by the table 
is characteristic of the damping behaviour of sandwich beans 
and plates with constrained damping layer treatment [8,24,91] . 

Table 3.2 gives the coincidence frequencies and loss 
factors for different values of the core loss factor p and 
for typical values of the core shear and geometric para- 
meters. It is observed that the loss factor varies linearly 
with respect to the core loss factor p, for the particular 
value of g considered. There is practically no change 
in the coincidence frequency with p. But marginal changes 
in the coincidence frequencies with p can be observed for 
greater values of g as will be discussed later in connec- 
tion with transmission loss curves. 

Table 3.3 gives the coincidence frequencies for dif- 
ferent values of the core geometric parameter Y, for the 
particular value of the shear parameter g* = 30.97(1+0.31). 

The coincidence frequencies increase with increase in Y 
for all values of 9 . However, the loss factor does not 
increase monotonically withY. Moreover, it is observed 
that the variation of the loss factor with 0 , for all 
values of Y considered is similar, in that the ratio of 
the loss factors for any two angles of incidence is very 
nearly a constant for all values of Y . 



Table 3.2 - Coincidence Frequencies and loss Factors for different values of p 
^ » 0.001 m, Y = 31.5 = 2,24 h, ) , g -JdFf 3o.^y- 
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Table 3.3 - Coincidence Frequencies and Loss Factors for Different values of Y 

b- = 0.001, i = p = 0.3 

•** 50.97 
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Figure 3.2 presents the variation of transmission loss 
with frequency for different values of the shear parameter 
and an incidence angle of e= 75°. The other parameters 
which are kept constant are shown in the figure. For all 
values of g the transmiss ion loss curves have dips at the 
corresponding coincidence frequencies. As had been already 
observed from Table 3.1, the coincidence frequencies decrease 
with increasing g which is also seen from the coincidence 
dips. For low values of g the transmission loss curves 
follow the mass law almost uptill the coincidence frequency 
and have sharp dips at the coincidence frequency. For high 
values of g too, the transmission loss curves follow the 
mass law at low frequencies, but start deviating from the 
mass law much earlier than the coincidence frequency. The 

less 

sound transmission's significantly higher at intermediate 
frequency ranges (1000 - 4000 Hz) for the values of g = 10 
and 100 as compared to those for g = 5000 and 10000. It 
is mainly due to the fact that the coincidence frequencies 
for the low values of g are much higher than for the lar- 
ger values of g and the transmission loss curves for the 
latter tend to droop early, while for the former they con- 
tinue to follow the mass law. But beyond the coincidence 
frequencies, the transmission loss increases rapidly with 
increase In frequency for the large values of the shear 
parameter g, while for the smaller values of g the 
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transmission loss decreases from the mass law as the cor- 
responding coincidence frequencies are approached. Hence 
at these frequencies the transmission loss is very much 
greater for large values of g, especially as the coin- 
cidence dips for small gs are sharp and deep resulting 
in a very low transmission loss at coincidence. The 
low values of coincidence transmission loss for low values 
of the shear parameter can be attributed to the s mal l 
values of the loss factors at coincidence corresponding 
to these shear parameters. 

It is also seen from the graphs presented that the 
coincidence transmission loss can be considerably improved 
(about 20 dB) by an optimal choice of the core shear para- 
meter g. For example, for the case of g = 1000, the trans- 
mission loss is the highest at coincidence, which is due 
to the fact that the coincidence loss factor for g = 1000 
is the highest for the angle of Incidence 0 = 75°» as can 
be seen from Table 3.1. The transmission loss curve for 
this value of g is very nearly flat and the coincidence 
dip is barely visible. Curves shown in Figure 3.2, thus 
clearly indicate that the transmission characteristics of 
the sandwich panel can be considerably improved, by a sui- 
table choice of the shear parameter g, depending on the 
frequency ranges of interest. 
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Figure 3.3 shows the transmission loss curves for 
various values of g for an incidence angle of © = 45°. 

The same observations are made as for e = 75°# except 
that the transmission loss values are higher and the coin- 
cidence frequencies are also higher. Another observation 
that can be made from the figure is that the optimum value 
of the shear parameter g for maximum transmission loss 
at coincidence for 9 = 45° need not be the same as for 
9 = 75°. 

Figure 3.4 shows the behaviour of the transmission 
loss curves with the core loss factor {3 , for two different 
values of the core shear parameter g. ( g = 30.97 and g = 5000 ). 
It is seen that for the low value of g, the effect of p on 
the transmission loss is minimal except at the coincidence 
frequency, where with increasing p , increasing transmission 
loss is obtained. But the coincidence dips are steep. The 
effect of p on the transmission loss is more discernible 
in the case g = 5000 . Significant increase in the coin- 
cidence transmission loss (13 to 17 dB) of the sandwich 
panel is obtained with increase in p. It is because, the 
effect of increase in p is compounded by the effect of 
increase in the core shear parameter resulting in high 
loss factors at coincidence as can be seen from Table 3.1. 
The coincidence dips are also spread out and flattened for 
the same reason. Another observation that is made from 
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the graphs is that , while increase in p has virtually no 
effect on the location of the coincidence frequencies for 
g = 30.97 f the coincidence frequencies are shifted to the 
left with increase in p for g = 5000. Beyond the coin- 
cidence frequency also, higher values of p have a beneficial 
effect on the transmission characteristics of the sandwich 
panel for the higher value of g. This is not so for g=30.97. 

The transmission loss curves for different values 
of the geometric parameter is shown in Figure 3.5 for a 

ft 

constant value of the shear modulus of the core (G- =2i .76x10° 

2 

B/m , p = 0.3). The geometric parameters correspond to 
core thicknesses of one half, equal and four times the thick- 
ness of the face plate respectively. As can be seen from 
the graphs, till coincidence, increase in core thickness 
results in increasing transmission loss. As have been 
already observed from Table 3.3 the effect of increase 
in the geometric parameter is to increase the coincidence 
frequency which can also be observed from the graphs. 

Figure 3.6 shows the variation of the transmission 
loss with incidence angle. As is obvious from equation 
(3.20), the transmission loss decreases with increasing 
angles of incidence and it is depicted in the figure. 

The coincidence frequency also decreases with incareasing 
angles of incidence which is also seen from the coin- 
cidence dips. 
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FIG. 3.5 VARIATION OF TL WITH GEOMETRIC PARAMETER 
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In Figure 3.7, the effect of increasing the face plate 
thickness Is shown. Simultaneously the core thickness is 
also increased in the same proportion, meaning thereby that 
the geometric parameter Y is kept constant. The shear 
parameter is also kept constant in these curves. With in- 
creasing thickness of the face plate the coincidence frequency 
is decreased. This is clear from equation (3.2 5). The 
coincidence frequency is proportional to the square root 
of the ratio of the surface mass density to the total bend- 
ing rigidity of the sandwich panel. While increases as 
the cube of the face plate thickness the surface mass den- 
sity p varies only linearly with face plate thickness. 

Hence the decrease in coincidence frequency with increasing 
h^. The transmission loss Increases with increase in h^, 
for in the mass law region the mass of the sandwich panel 
is considerably increased. The coincidence transmission 
loss is also higher for greater thickness of the face plate. 
This may be due to increased loss factors at coincidence 
with increase in face plate thickness. 

Figure 3.8 compares the transmission loss of an equi- 
valent homogeneous plate having the same surface mass den- 
sity of the sandwich panel and of the face plate material 
with the transmission loss of the sandwich panel for two 
typical values of the core shear parameter g and at an 
incidence angle of 75°. It is evident from the graphs 
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that the sandwich helps in shifting the coincidence fre- 
quencies towards higher frequencies. Also, because of the 
damping layer the transmission loss at coincidence is in- 
creased to significant levels. 

3.5 COHCHJSIOMS 

The problem of the sound transmission loss through 
an infinite sandwich panel with a constrained damping layer 
treatment is formulated in this chapter. Expressions for 
the sound trams miss ion loss and the coincidence frequency 
are derived in terms of the sandwich core parameters and the 
incidence angle. From a parametric study of the transmission 
loss of the sandwich panel by the following conclusions can 
be drawn. 

i) The transmission loss is more sensitive to the variation 
of the core shear parameter g than any other parameter. 

ii) High coincidence transmission loss is obtained by a sui- 
table choice of g and a high value of the core loss 
factor p . 

iii) The transmission loss of the sandwich panel increases 
with the thickness of the core below the coincidence 
frequency. 

iv) The coincidence frequency reduces significantly with 
increase in g , while it increases with increase 



in the geometric parameter Y. 

) The sandwich panel has improved sound transmis 
sion characteristics as compared to a homogene 
ous panel of equivalent surface mass density. 



CHAPTER 4 


NOISE TRANSMISSION AND STRUCTURAL RESPONSE OF DAMPED FINITE 
SANDWICH PANELS BACKED BY RECTANGULAR CAVITY 

4.1 INTRODUCTION 

Several works have been presented in the past concern- 
ing the behaviour of cavity backed finite panels subjected 
to external airborne excitation. The evaluations of Pret— 
love [49], Dowell [50], Bhattacharya and Crocker [52], Guy 
and Bhattacharya [53], and Guy [56] have led to a progres- 
sively greater level of understanding of the phenomenon and 
improved accuracy of prediction. Dowell, et al [55] in a 
recent paper had considered the interaction between the 
internal sound pressure field and the elastic flexible wall 
of an enclosure in an 1 a cous to elasticity * formulation. They 
adopted a series expansion in terms of the structural normal 
modes and the cavity acoustic modes and reduced the problem 
of determining the structural response and the sound pres- 
sure inside the enclosure to a matrix eigenvalue problem with 
gyroscoping coupling. All the works mentioned above dealt 
with elastic homogeneous panels and assumed simply supported 
panel edge conditions. Moreover, the analysis was mainly 
restricted to external harmonic excitation. These investi- 
gations led to interesting findings with respect to resonance 
location of the plate-cavity system and detailing the 
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mechanism of energy transfer. The problem was also investi- 
gated by use of statistical energy analysis by Lyon [57] and 
Eichler [58] in the high frequency ranges. 

Recently, Vaicaitis [2] considered the noise trans- 
mission through a sandwich plate with a constrained visco- 
elastic damping layer treatment into a closed rectangular 
cavity, for a stationary random external excitation. He 
used a normal mode analysis, and showed that significant 
noise reduction can be achieved in the frequency range of 
0 - 1000 Hz, by the additive damping treatment. In his 
analysis also the edges of the plate were assumed to be 
simply supported. Another assumption male in his analysis 
was that the plate motion was unaffected by the cavity pres- 
sure and could be determined by the in-vacuo flexural res- 
ponse of the plate. But the effect of the plate motion on 
the cavity pressure was taken into consideration. Thus 
only a partial coupling between the structure and the 
acoustic field was considered in his analysis. 

This chapter concerns with the problem of noise trans- 
mission and structural response of elastically supported 
sandwich panels with constrained damping layer treatment 
backed by a cavity. Both harmonic and stationary external 
random pressure excitations are considered. Initially, an 
analysis based on Vaicaitis [2] is presented, where only 
a partial coupling of the structural motion and internal 
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cavity pressure is taken into account. Later, a more 
general analysis, allowing a full interaction between the 
internal sound field and the structural motion is made in 
an acoustoelastic formulation. 

The y~wise edges at x = + a/2 (Figure 4.1) are assumed 
to be elastically supported and the x-wise edges are assu- 
med to be simply supported. These boundary conditions re- 
present realistic boundary conditions that would obtain in 
aircraft skin-stringer constructions. The justification for 
the above assumption have been explained earlier in Chapter 2. 

Interior cabin noise in propeller driven aircrafts, has 
maximum intensities in the low frequency range (0 - 1000 Hz) 
corresponding to the blade passing frequencies. In this 
analysis the noise transmission studies are restricted to 
this frequency range . 

The response and noise transmission calculations in 
both the formulations use the forced damped normal mode 
analysis in view of the arbitrary boundary conditions at 
the panel edges and arbitrary damping. 

In the acoustoelastic formulation the problem can be 
reduced to a matrix eigenvalue problem with gyroscopic 
coupling. Determining the eigenvalues and the eigenvectors 
for the combined natural frequencies of the cavity -panel 
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system involves cumbersome complex arithmetic due to the 
presence of the damping layer, whose shear modulus is 
represented by a complex number. However, since the con- 
cern here is mainly with response and noise transmission 
characteristics, the eigenvalue problem is not solved di- 
rectly. The response and noise transmission are obtained 
by a direct matrix inversion scheme, taking advantage of 
the diagonal nature of certain matrices in a partition ma- 
trix analysis. 

Parametric studies are made with respect to the sand- 
wich core parameters and computer investigations show that 
the damping layer improves the noise reduction correspond- 
ing to structural resonances, especially near the first 
fundamental natural frequency of the panel. The damping 
layer significantly reduces the structural response levels 
near the structural resonant frequencies. 

4.2 FORMULATION OF THE NOISE TRANSMISSION PROBLM NEGLECTING 
THE EFFECT OF CAVITY ON STRUCTXRAL MOTION 

4.2.1 Acoustic Solution 

Consider the rectangular enclosure shown in Figure 4.2. 
The wall at z = O is flexible being the damped sandwich panel 
of concern. The other walls of the enclosure are assumed to 
be acoustically rigid. 
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FIG. 4.2 RECTANGULAR CAVITY - SANDWICH PANEL 



82 


The acoustic pressure p c inside the cavity is governed 
by the wave equation 



p c 



with 


ap a 2 w 

~ — = - P — 2 on K 

a n at ? 


(4.1) 


(4.2) 


=0 on Ag 

2 a 2 q 2 

where, = r- + — — n + — « , p and c are the 

d x 2 dz 2 

^3-U-i.d density and the local acoustic velocity in the medium 
respectively. w is the displacement of the flexible wall 
in the normal direction n (z direction). Ap and A R indicate 
the flexible and rigid portions of the enclosure respecti- 
vely and t represents time. 


It is assumed in this particular analysis, that the 
cavity pressure does not affect the wall motion and the 
radiation loading on the plate surface, external to the 
cavity is negligible. 

The solution of equation (4.1) can be written as 

oo oo 

P G (x,y,z,t ) =2 2 C (z,t) f .(x,y) (4.3) 

° i*0 3=0 13 ±3 

where, f^ are the cavity modes with rigid walls at 
x = + a/2 and y = 0,b (Figure 4.2) given by f^=cos^-(x+ j|)cos^~2,. 



'83 


The boundary conditions for the pressure coefficients 
C. . in equation (4.3) are non-homogeneous. The solution 

i tj 

for Ch .. can be effected by transformation of the homogeneous 
differential equation with non-homogeneous boundary condi- 
tions into a non-homogeneous differential equation with 
homogeneous boundary conditions [92]. 

Then one can express, 

°i3 = + Z ij <*•*> 

where, a. - is the solution of the associated homogeneous 

XJ 

problem and Z. . is chosen to satisfy the given boundary 

X J 

conditions . 


Substitution of equation (4.4) into equation (4.3) and 
(4.1) gives 


a 2 a ± . • 1 1 

- J “ij - A ij “ij = V -13 


s 2 z. 


Z_. 1 - 


u + 


dz 




(4.5) 


where the number of dots represents the order of differen- 

2 2 2 2 2 

tiation with respect to t and A . • = % (i /a + 3 /b ). 

XJ 


Since Z. . is the solution satisfying the boundary 
13 


condit ions , 
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az. 


41 


3z 


2=0 




(4.6) 


co oo 


where, 2 z G .(t) f .(x,y) = - P» (x,y,0,t) 

i=0 j=0 3 0 


Using the orthogonality property of the cavity modes 


fl 3 U,y>, it can he shown that 


e. . 


a /2 b 


G. i (t) f / p w f . (x,y ) dxdy 

1J ab -a /2 o 13 


where , 


e. .= 
13 


1 for i = 0 and 3=0 

2 i / 0 or 3 / 0 

4 i / 0 and j / 0 


(4.7) 


can now be chosen as, 
Z ± . = (z - z 2 / 2 d )G ± j 


(4.8) 


where, d is the cavity depth. It can be seen that the form 
of equation (4.8) satisfies the boundary conditions of equa- 
tion ( 4 . 6 ). 


From equation (4.5)* the solution for a^.( 2 ,t ) can 
be obtained as. 


oo 

a ±;j (z,t) = z B iik (t) cos 


krcz 


k=0 


ijk' 


(4.9) 
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where the B-^’s are governed by a single degree of freedom 
system equation of the form 


B i jk + 2 e ijk “ ijk B ijk + 40 ijk B ijk ~ B ijk 


( 4 . 10 ) 


where is an equivalent viscous modal damping coefficient 

which accounts for the absorption losses in the wall and 
air damping in the cavity. The modal damping coefficient in 
the analysis is taken to be. 


£ — £ ( Oft / (jJ V 

e i;jk ” ool ^ ool 7 ijk' 


(4.11) 


based on the works of references [2,5]. 

B ijk in e l ua ' t i° n (4.10) is the generalized force correspond- 

-i- T_ 

ing to the ijk 1 cavity mode, given by 

t — 

ijk “ 


2 2 
- c e, d 1 .. 3 Z . . k-ftz 

P-,.^ = - f( — £ Z. -+1^ - Z. ^ - f)oos 

d. o c j i j d z*~ d 


dz 


(4.12) 


where, e fc 


1 for k = 0 

2 for k ^ 0 


and w. are the s quare of the cavity resonant frequencies 
given by 


.2 2 r , 2 2 / , 2 \ 

“ijk = 0 (A u + k % /& ) 


(4.13) 


Prom equations (4.9), (4.4) and (4.3), the cavity pressure 
is given by 
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knz 

d 

(4.14) 

4.2.2 Structural Response 

&) Harmonic loading ; Consider a plane harmonic pressure 
wave incident on the plate at an angle 0 to the normal. The 
external pressure field can he expressed as 

P e U,y,Q,t) = p Q exp [i(« t-k’xsin e)3 (4.15) 

where k' =ti) /c is the acoustic wave number. The pressure 
distribution is assumed to be uniform along the y direction. 
As it has been assumed in the analysis that the plate motion 
is not affected by the cavity pressure and the radiation 
loading on the plate external to the cavity is negligible, 
the structural response can be directly obtained as the 
in— vacuo structural motion. A forced damped normal mode 
analysis outlined in Chapter 2 is used for the same. 

OO OG 

Expressing w(x,y,t) = z Z W _ n (x,y) a (t) 

m=l n=l 

( 4 . 16 ) 

where W mn (x,y) are f° rce d damped normal modes, the 

damped normal coordinates are given by 


OO' 'OO OO 

p (x,y,z,t) = Z Z Z B (t) f..(x,y ) cos 
w i=0 j=0 k=0 


OO OO 


+ . 2 2 Z ±i (z,t ) f (x,y) 

1=0 j =0 ^ 
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W*> 


P ( t ) 

r -nn v ' 

Tan 


(4.17) 


where, p Qn and p rar> are given by equat ions (2.31 ) and (2.32 ) 
respectively . 


The particular forms of p Tr>n and for the different 

n2 

boundary conditions are given in Appendix B. (d ( 1+in n ) 

4-TU 

is the complex natural frequency in the mn damped mode 
and t} ^ is the corresponding modal loss factor. 


b) Random loading ; Consider a random pressure excitation 
on the plate. The pressure field is assumed to be stationary 
with respect to time and homogeneous with respect to the 
spatial indexing parameters. The power spectral density of 
the external pressure can be expressed as in equation (2.33). 

S pp (3 l’ y l ; 3 2' y 2 i “ ) = S P 0 P 0 ( “ )r xy (x 2~*l ; y 2- y li“> > 

(4.18) 

In further analysis, the external excitation is 

assumed to be white noise with S ( u> ) = S , and uniform 

p o p o ° 

with respect to the spatial co-ordinates , that is, 
r xy^ x L ,y l* x 2»y2* w ^ = 1#0 * Tiien tiie power spectral density 
of the response is given by 

OO OO CO OO 

^ww^ x l ,y l’ 2 2 ,y 2* ^ ^ >y l ^ ^rs^ x 2 ,y 2^ * 

x m=l n=l r=sl s=l 


H ( w)H* (w) S I f <d) 
mn v ' rs' 3 o mars' 


(4.19) 
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The various quantities in equation (4.19) are defined 
by equation (2.35) and (2.36) and for the different plate 
boundary conditions are given in Appendix B. 


4 .2.3 Noise Transmission into Cavity 

a) Harmonic loading ; Having determined the structural res- 
ponse, the pressure inside the cavity can be calculated in 
the following way. Substituting equation (4.16) into equa- 

4* 

tion (4.8), one gets. 




p e . . oo co .. 

“* £ £ 1 . . q_ 

ab m=l nil ^ 


(4.20) 


where , 

a/2 b 

Wr-/ ! \n (x ’ y) ay4x (4 - 21) 

—a/2 o 


Prom equat ions ( 4 . 9 ), (4.11), (4.12) and (4.20) it 
follows that. 


2 

pc d e. . 

Bj = - 3J. 6j 


OO oo 


ijk 
where , 


ab 


6 ijk = 


(■*£ - A ij )/3 ' 1/d ' 


« k mil nil ^ 
2 






w ijk“ w+21e ijk w iak u 


(A.. - a) 2 / c 2 ) 


-] for* k = 0 


~ k ^ 2 (w ijk- W +2i£ iok "ijk^ 


for k 4 - 0 


(4.22) 


(4.23) 
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The pressure inside the cavity can new be obtained 
from equation (4.14) as 


P c (x,y,z,t ) 



s 

i=0 


CC CO 

t y 
j=0 k=0 


e ij 6 ijk 


E. 

m=l 


OO 

2 

n=l 


^knij %n' 


f 1;j (x,y) cos - ( z -z 2 /2d) X. 


OO OO 

E E 


OO 03 


i=G j=0 i 3 


2 2 L . . a f. . (x,y) 

m=l n=l mni ^ ^ 


The noise reduction (HR) through the sandwich panel 
can be calculated by computing the cavity pressure at the 
panel at z = 0. The noise reduction through the centre of 
the panel (x = 0, y = b/2, z = 0) is then expressed as. 


HR = 20 log ?2_ (4.25) 

jP c (°> 2 * °)] 

b) Random loading : Prom equations (4.8) and (4.19) the 
cross power spectral density between the functions G-^ and 
G lo is given by 

2 4 

• p e. .e, a) 00 °° 00 00 

((0) = S r E E E E 

^ijlo a d 0 m=l n=l r=l s=l 

H mn Wj H rs 4lo ] < 4 ‘ 26) 

Prom equations (4.14) and (4.22) the cross spectral 
density of the cavity pressures at (x = y = y^, z = 0) 
and at (x = Xg, y = y 2 » 2 = 0) can be expressed as. 
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S , . ,, , 4,2 r £ Z Z 2 2 Z 

p c p c^ x l ,y l’ x 2 ,y 2 # ^~ c d [± j k 1 o h 


k 6 loh f i 3 < * y :i. 5 Ro ( *2 > y 2 )s e ljlo ( ] 

(4.27) 

The auto power spectral density of the cavity pressure 
at a point (x,y ,0) in the sandwich panel can he obtained 
from equation (4.27) by setting 2 ^ = Xg = x and y^ = y 2 =y . 

The noise reduction at the centre of the plate is 
calculated according to the formula 


NR = 10 log 


S (0, b/2, 0,(i) ) 

^(rC 


(4.28) 


4.3 ACQUSTOELASTIC FORMULATION 
4.3.1 Acoustic solution 

Consider again the cavity shown in Fig. (4. 2) occupying 
a volume V = abd. As before, it is assumed that the wall 
formed by the damped sandwich panel at z = 0 is flexible 
and the other walls are acoustically rigid. 


The cavity pressure p is governed by the following 

C 


e quat ions 


_,2 

V P, 


L2< 

av 


0 


(4.1) 
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0P C 0 2 W 

= -p— on 

3n dir 

= 0 on 

Equation (4.1) for 


■% 

homogeneous 


8 Pc 

= 0, has normal mode solutions 

8 n 

the equations 

V 2l? n = -< V /c)2 


(4.2) 

boundary conditions 
3? n eiWt which satisfy 

(4.29) 


= 0 on A (4.30) 

an 

where, is the n^ 11 acoustic natural frequency, P n is 

the corresponding cavity mode, A is the total surface area 
of the enclosure. Rote that the index n contains the 
indices i,j,k of the previous section in a compressed form 
which represent the modal shapes in the x,y and z directions. 
(3? n corresponds to x,y p) of the previous section). 

The orthogonality condition of the acoustic normal 
modes can be expressed as. 



f 


F„ F_ dV = M. 6_„ 
r n nA rn 


where 6 rn represents the Kronecker delta. 


(4.31) 


The wave equation (4.1) can be transformed into a 
set of ordinary differential equations by using Green’s 
theorem in the form. 
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t. o o r 3P n tj, dp 

/ (P V F n ~ F n V P) dV = / (p P n — )dA 

Y n n A 3n 9n 

(4.32) 


The acoustic pressure inside the cavity can be expan- 
ded in terms of the acoustic normal modes P as 

n 




P P 
n n 


M, 


nA 


(4.33) 


The pressure coefficients P n in the acoustic normal 
mode expansion are governed by the set of ordinary differ- 
ential equations 


P + w , p _ _ 

n nA n 



(4.34) 


where. 



1 

“T~: 

pc Y 


f 

Y 


p P dV 
*c n 


(4.35) 


and 



1 

/ 

A-p Ap 


w F n iA P 


(4.36) 


Wall absorption and air damping in the cavity can be expres- 
sed in terms of an equivalent viscous modal damping in 
which case equation (4.34) can be written in the form 


Ap »• 

P + 2e W p + w 2 F n = “ Y~ W n 
r n + ^ e nA nA n nA 


(4.37) 


where, is the modal damping coefficient in the n 


th 


acoustic mode. 
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4.3.2 Coupled Acoustoelastic Equations 

The equation of motion for the flexible sandwich panel 
with a constrained viscoelastic damping layer is given by 
equation (2.9) 

2 

y 6 w-g* (1+Y ) V 4 w + ~ £-2 [ V 2 w-g'w] = — ( \7 2 p-S'p) 

"t at D t 

(4.38) 

In the above equation p is the net pressure 

excitation on the plate -given by p = p - p , where p 

is the internal cavity pressure and p is the external 

pressure loading on the plate. The radiation loading on 

the plate surface external to the cavity is neglected. 

* 

As before, the response of the damped sandwich plate 
for a harmonic or stationary random pressure excitation can 
be expanded in terms of forced damped normal modes in the 
form, 

' = ( 4 - 59 ) 

m 

Here, in the index m are contained the indices of 
both m,n corresponding to the mode shapes in the x and 
y directions. In the terminology of equation (4.16), 
corresponds to q^ n and <p m the damped normal mode corres- 
ponds to W mri (x,y ). Using the method of damped normal 
mode analysis outlined in Chapter 2, the damped normal 
co-ordinates q^s are governed by the ordinary differential 
equations 
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•* yjO 

+ %s = Snc + Sae (4 ' 40) 

where M^g is the generalized mass corresponding to the 
damped normal mode given by , 

M ms = f ^ <4 dA (4.41) 

% 

%c = / Pc dA and Sae = / p e 9 m dA ^ the Senerali- 

A F 

zed forces corresponding to the cavity and external pres- 
sures respectively. It should be noted that in general, 

M ms* Snc and Sue are com P lez values. Using equations 
(4.33) to (4.36) and equations (4.39) to (4.41) the coupled 
a cous to— structural equations can be expressed as 

+ 2e nA “nA K + “L. P n = ~ 2 ^ V < 4 - 42 > 

m v 

M ms^*^n + 0) ms^ 1+il Tms^9 m ] = p c •% s (4.43) 

n nA 

where, 

L nm = Aj F n ^m dA (4.44) 

Note that, L m corresponds to 1 ^. ^ of equation (4.21) 
divided by the area of the panel. 

4.3*3 Solution of Ac ous to elastic Problem 

It is shown by Dowell [5], that expressing equations 
(4.42) and (4.43) in terms of the acoustic potential and 
in matrix notation, and limiting the summation in the 
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equations to a finite number of acoustic and structural 
normal modes, the problem reduces to a matrix eigenvalue 
problem with gyroscopic coupling, as would arise in the 
dynamics of spinning structural systems. The solution of 
the eigenvalue problem will give the coupled natural 
frequencies of the struct tare-cavity system. The eigen- 
vectors can be used to obtain the pressure coefficients 
3?n s and the structural coordinates q ffl s . 

The solution of the eigenvalue problem for gyroscopic 
systems can be obtained by expressing the coupled acousto- 
elastic equations in state space form, in terms of displace- 
ments and velocity coordinates. Then one has to deal with 
2 (ntf-n) first order ordinary differential equations defined 
by two matrices one symmetric and the other skew symmetric. 
Meirovit ch [93] has proposed a new method of solution of 
such problems by taking advantage of the special nature 
of the problem. But his method is applicable only if the 
elements in the symmetric and skew symmetric matrices are 
real. 


However, in the present work, many of the matrix 
elements are complex valued, due to the presence of the 
additive damping layer. Hence an eigenvalue analysis 
of the problem is a difficult proposition, because of 
the cumbersome complex arithmetic and the size of the 



matrices involved. Moreover, the main interest in this 
analyst's is to obtain the noise reduction inside the ca- 
vity and the structural response . They can be determined 

d ire ctly by a matrix inversion procedure, which takes ad- 

• 

vantage of the diagonal nature of some of the matrices in 
the analysis. Equations (4.42) and (4.43) can be expres- 
sed in matrix form as. 


"t 1 ! ! ' 

(nXn) , (M.) 

] 

J 

(& ~ 
(nxl) 


[- 2e . w ■>]» 

L nA nA J j 

(nxn) , 

[0] ‘ 
(nxm) 


1*1 | 
(nxl) 

J 

A 

— — 

y + 



— — . 

1 

> 

to] ; [i] 


Jq} 

f ' 

[ 

[0] 

[0] 


■£} f 

(mxn) , (mxm) 


Cmxl) 

J 


(mxn) * 

L ! 

(mxm) 


(mxl V 


! to] 

’j 

;• W ' 


p 0 ^ 

(nxn) j (nxm) 


(nxl) 


t B 2 ] l C" , ^s( 1 +lt Ws ).] 

(mxn) , (mxm) 

| 

{ 4 ! | 

> = 

w f 


' ( mxl )j 

f 

(_(mxl ) 


where, [I] and [0] are identity and null matrices of appro- 
priate orders indicated. Typical elements of the matrices 
[B^] and [Bg] are given by. 


F 


li j 


13 


and 


3 2ij 


p c A-p 


M - , M- „ 
3A is 


3i 


and ^q| are the vectors of the pressure coefficients 
P n s and generalized coordinates and |Q|is the external 


force vector whose m element is given by 


ms 
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The letters Inside the brackets indicate the order of the 

* 

mat rice- n being the number of acoustic cavity modes, m, 
the number of structural damped normal modes considered in 
the modal summations. 


a) Harmonic loading : Consider an external pressure excita- 
tion with a harmonic component e* ^ Matrix equations(4.45 ) 
are similar to equations of motion of linear multi-degree 
of freedom vibrating systems with both inertial and elas- 
tic coupling. The steady state solutions of such a system 
are harmoni c with the same frequency to . Hence assuming 

{E} = {p} ana and = {q} e iu * 

S' 

e quat ion (4.4#) can be written as 



where, the matrix [B] is given by. 


(4.46 ) 



_ _ _ _ _ 

t%] ! 


(4.47) 


In equation ( 4 . 47 ) [^ Z n ^( w ) -] and [" Z^ s ( to ) -] are diagonal 
matrices having diagonal elements 


Z nA<“> =“L + 21 E nA “aA - “ 2 

=“ms ' 1+ lr W> 


(4.48) 

(4.49) 
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Z nA (W) ’ and 2^ s ( w ) can be respectively regarded as the 
complex cavity impedance function in the n a cavity mode 
and the structural impedance function in the m damped 
normal mode. 


The complex amplitudes |p| and a can be obtained 


as , 



(4.50) 


The pressure inside the cavity can be obtained from 
equation (4.33) as. 


p 0 = P o 2 e^ S VAl 

n 


(4.51) 


and the structural response at a point (x,y) on the sur- 
face of the plate as, 


/ \ jjGOt 

w(x,y) = e S 71 tp_(x,y ) 
mm 


(4.52) 


The noise reduction inside the cavity from equation 
(4.28) is 


NR = 20 log 


fc2 f f 


(4.53) 


where, p Q is the amplitude of the external harmonic pres- 


sure excitation. 
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b. Random loading : Consider the external excitation to 
be a stationary random process with power spectral den- 
sity function, given by equation (4.18). The cross spec- 
tral density between the generalized forces and is 
then given by. 




S ( a)) j (co ) 

PqPq 

M ms M rs 


(4.54) 


where , 




J / r ^ x l ,y l» x 2 ,y 2 ;a5 ) ! Pm^ X L ,y 1 "> dA 2. d - A 2 

"p % 

(4.55) 


Note that ) in the foregoing equation corresponds to 

the acceptance functions I^g defined in equation (4.19), 
the index m including both the indices m and n and the 
index r including both the indices r and s in a compres- 
sed form. The asterisk indicates complex conjugation. 


The cross spectral density matrix of the generalized 
forces can be constructed from equations (4.54) and (4.55) 
as 


[s QQ (®)] = S p 0 P 0 (a)) C I(03) ] (4 * 56) 

where a typical element of the acceptance matrix [I] is 
given by 


I „( w) 
I mr (o>) 

M ms M rs 


(4.57) 
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The cross spectral density matrix of the pressure 
coefficients s and the structural generalized coor- 
dinates q^s can be shown to be 


[S^tw)] = 

(iH-m)x(nHn) 


[Spp(“)] 

(nxn) 


i - 


~i 


[S 


qp( u )l 

(mxn) 


CS ?q «o)] 

(nxn ) 


[ S qo ( “ )] 

(mm) 


[A] 


[0] I [0] 


CO] ; [SqqC-H 


[a*] 


T 


(4.58) 


where the steers cript T refers to transpose and the 
asterisk again represents complex conjugation , 


Prom equations (4.33) and (4.58) the power spectral 
density of the cavity pressure at any point in the cavity 
is. 


S P C P =P 2 c 4 5 [S p p(u))] ja*| (4.59) 

where, is a n-dimens ional vector whose n element 

is F n U,y,z)/M nA . 


Similarly using equations (4.16) and (4.58) the 
power spectral density of the displacement w of the 
panel at a point (x,y) is given by 


S w (x,y,o.) = [S qq (»)] 


(4.60) 


th 


where, |y| i s a m-d imensional vector whose m element is 

<p m (x,y). 
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The noise transmission characteristic of the sandwich 
panel can he expressed in terms of the noise reduction 


m = 10 log [s poPo ( x , y , 0 , w)/ ^ (*,7.0.-)] 


( 4 . 61 ) 


4 * 3 * 4 ISZersion of fB] K&trix 


Thus both for harmonic and random external acoustic 
excitation, the structural response, the cavity pressure 
and the noise transmission characteristics of the sandwich 
panel can be obtained by the inversion of the matrix [B], 


It is not necessary to invert the whole matrix [B] 
which is of order (in-m) x (n+m). Advantage can be taken of 
the presence of the diagonal matrices and 

£ ^ns ( w ) -1 • Writing [B] 1 = [a] in partitioned form. 


[A] = 


[A X ] 

(nxn) 

C ^3 

(mxn) 

it can be shown that. 


t - 


[As] 

(nxn) 

[a 4 ] 

(mxm) J 


( 4 . 62 ) 


-“ 2 [ a ,_] [' Z ^-]" 1 [%] 
C V = ~" 2 [%] c-s^- 1 [Bl ] 

[Ag] = ry [a 4 ] 

CBg] [AJ 


( 4 . 63 ) 
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[~ Z ^x -.] and [~ Z^l] are just the reciprocals of the 
diagonal net rices [- Z n A -] and [~ 2^ -] respectively. 

Hence, the problem of determining the inverse of matrix 
[B] of order (rn-m) x (n+n) reduces to finding the inver- 
ses of two matrices of order (n x n) and (m x m) respec- 
t ively . 

It is further shown here, that it is sufficient to find 
only the matrix [A^] which is the inverse of a matrix of 
order (mxm), for the determination of the cavity pressure 
and the structural response and not the invers es of two 
matrices as indicated above. It is so because, from equa- 
tion (4.58), it is seen that the spectral density matrices 
are given by, 

[Spp(u) ] = [Aj] [S QQ ] [Ag*] 1 

[S 49 («> 1 - u 4 ] [S QQ ] (4 _ 64) 

tV u> 1 ■ [a 2 ] cO 1 . 

[S qp (0.) ] = [A 4 ] [S Qq ] [Ap*] 1 

Only the matrices [Ag] and [A.] are needed for a 
complete solution of the problem. Prom equation (4.65) 

[Ag] is obtained from [A^] by a simple transformation 
and so the solution of the acoustoelastic problem re- 
quires the inversion of the matrix of order (mxm). 
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4.4 APHtOZH-IATB F 03MULA1 IOHS 

Before the discussion of the results obtained for a 
typical cavity-sandwich panel system, it is purposeful to 
consider certain simplified models as considered by Dowell[5] , 
that would lead to a better understanding of the results. 


If the exciting frequency w , say w g , is very nearly 
equal to a structural resonant frequency 0 )^ , the coupled 
acousto elastic equation in the structural coordinates can be 
decoupled by assuming that the main contribution to the struc- 
tural response will arise due to q^ alone. Hence from equa- 
tions (4.42) for harmonic loading 


Ar 


P n = 


nm 


Sn 


(— u) + 0 ) .+2ie .oj , a) * ) 
ms nA nA nA ms 


(4.65) 


and from equation (4.43), 


P c 2 A 2 

^ms^^m ms^ + ^ms = Z 


[ ^ 


.2 

J nm 


L - V 


11 M nA<“ l!S + “nA +2ie nA “nA “mV 

3 B (4.66) 


Prom equation (4.66), it is seen that cavity modes 
having frequencies greater than w mg add to the inertia 
of the panel and cavity modes having frequencies less than 
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m add to the stiffness of the plate. The n = 0 mode 
ms 

6) = 0 always adds to the stiffness of the plate. 

%/■& 


If the structural resonant frequency « ^ and the 
cavity resonant frequencies are well separated, then the 
structural res onant frequencies will not be influenced by 
the cavity. Then in equation (4.66), the sum on the right 
hand side can be neglected as small compared to the external 
excitation term. Then the response to harmonic excitation 
will be governed by the modal loss factor of the structure 

i w-t 




i M_„rj 


xns^ms ms 


(4.67) 


Prcsa equations (4.33), (4.65) and (4.67), the pressure 
inside the cavity may be computed as. 


P c = 


2 . 

P c A p 

iT W v 


s 

n 


P 1 
n nm 

~W 


SmE 


M _M ,(- + u) 2 +2ie w . w* ) 

ms nA ms nA nA nA ms' 


(4.68) 


*2 

If w „ is very much less than the lowest nonzero acoustic 
ms 

resonance, only the n = 0 acoustic mode will be dominant and 
the cavity pressure p can be approximated as 


P c \ g o L om SaE 
V 1 ! M oA M ms 


(4.69) 


Thus from equation (4.69), under the conditions mentioned 
above, the cavity pressure is inversely proportional to 
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rj mS which can be likened to the energy dissipation 
parameter described by Mead [24]. The cavity pressure 
can at times exceed the external pressure giving rise to 

9t-p 

negative noise reduction for low values of ^ mS . It 
can also be observed that the noise transmission charac- 
teristics of the panel at structural resonances can be 
improved by increasing modal loss factors, especially for 
the fundamental structural mode and other low order modes. 


Significant changes in structural resonant frequen- 
cies will occur, mainly if a structural mode couples effec- 
tively with the n = 0 acoustic cavity mode. Equation (4.66) 
can then be approximated by 


M 


>>2 


ms 


[*^m + ^ms ^ 1+iri ms^~ ~ 


pc?4 £ % ♦ * 


i v 


om 


(4.70) 


From equation (4.70) it can be seen that the coupled 
structure-cavity frequency can be expressed as 


2 

u cs 


03 




ms 


4 * 


2 2 2 
p° % Sa 

v 7 


(4.71) 


Another case where the structural frequency may be 
changed due to the effects of cavity is when a structural 
resonant frequency and a cavity modal frequency coincide. 
In this case, there is a strong coupling between the struc- 
ture and the cavity and large sound pressure levels in the 
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cavity and large structural response corresponding to the 
coupled natural frequency of the cavity-structure system 
result. The phenomenon is described as multiple resonan- 
ces or double resonance in literature [53,55,90]. The re- 
sonant condition arises mainly because of an impedance 
matching between the cavity and structural impedances. Large 
press vires inside the cavity will also arise when two cavity 
frequencies are in close proximity . 

4.5 RESULTS AM? DISCUSSION' 

The following geometrical and material properties for 

the sandwich panel-cavity system are adopted in the computer 

investigations. The parameters which remain constant in 

most of the computations are. 

Thickness of face plate 

Length of panel along x direction 

Length of panel along y direction 

Youngs modulus of the face plate 
material E (aluminium) 

Density of face plate material 

h 

( aluminium ) 

Local velocity of sound in air 
Density of air 

Density of viscoelastic damping 
layer 

■ % 


= 0.0005 m 
a = 0.25 m 
b = 0. 5 m 


= 7.2 x 10 10 N/m 2 

p =2770 kgm/m 5 
3? 


c = 330 m/sec. 
p = 1.225 kgm/m 5 


p c = 0.2 p 
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For most of the studies, the depth of the cavity is 
taken to be d = 0.762 m. However, the effect of the 
cavity dept h on the structural response and noise reduc- 
tion inside the cavity are investigated by varying the 
cavity depths. As remarked earlier, in all cases the x-wise 
edges at y = 0,b are assumed to be simply supported (Figure 4.1). 

4.5.1 Order of Modes in Summation 

The boundary c end it ions at x = + a/2, are assumed to be 
identical and elastically supported. The external pressure 
excitation in the case of harmonic loading is assumed to 
be uniform. In the case of random loading the pressure 
excitation is assumed to be white noise with respect to the 
temporal parameter and spatially uniform with r (x^,y^; 
x 2» y 2’ W ^ = 1 *°* Under these assumptions, the contribu- 
tion to the generalized force in the case of harmonic load- 
ing or to the power spectrum of the generalized force in 
the case of random loading will arise only due to the sym- 
metric structural modes. Hence only the symmetric structural 
modes are considered in the summation. From equation (B.5) 
of Appendix B it is seen that the structural acoustic coup- 
ling terms represented by L mr ^ j have nonzero values only 
for the combinations of (m+i) and (n+j) being odd integers, 
for the simply sup ported boundary conditions. The following 
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combinations, (i = even, j = even, m = odd, n = odd), 

(i = even, j = odd, m = odd, n = even), (i = odd, j = even, 
n = even, n = odd), (i = odd, j = odd, m = even, n = even) 
are then admissible. But from the earlier observation, that 
the generalized forces have nonzero values only for the sym- 
metrical structural modes , the response and noise transmis- 
sion are computed by considering only the combinations 
(i = even, j = even, n = odd (symmetric), n = odd) of 
acoustic and structural modes. This entails no error in the 
formulation where the effect of the cavity on the structural 
motion is neglected. In the acoustoelastic formulation, it 
amounts to neglecting the contributions from the coupling 
between the odd order cavity modes and the even order struc- 
tural modes. The contributions from these tern© will any 
way be negligible compared to the terms considered in the 
summation because of the vanishing of the generalized forces 
corresponding to the even order structural modes represented 
on the right side of the equation (4.45). Since the fre- 
quency range of interest in this study has been restricted 
to 1000 Hz, the modal summations and the acoustic cavity 
modes include admissible combinations in them upto a fre- 
quency of 1000 Hz. Nine structural modes (m=l,3,5, n=l,3,5) 
and sixteen cavity modes (i=0,2, j =0,2, k=0, 1,2,3 ) are inclu- 
ded in the results presented. 
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4.5.2 Resonant Frequencies and loss Factors . 

The natural frequencies and the corresponding loss 
factors in the symmetrical modes, for different elastic 
support conditions at the edges x * + a/2 are presented 
in Table 4*1. It can be observed that the lowest natural 
frequency for transversely restrained edges could be less 
than the fundamental frequency for simply supported edges 
at x = + a/2, depending on the relative value of the strin- 
ger flexural stiffness parameter It has been noticed 

during the computer investigations, that increasing the 
value of to the order of 10^, results in natural fre- 
quencies add loss factors for the transversely restrained 
edges approaching the values of fixed edges. The natural 
frequencies and loss factors for the rotational ly restrai- 
ned edges with and without rivets do not differ very much 
for the particular values of the torsional stiffness para- 
meters K and K chosen and for the low value of the shear 

Sv 

parameter g* = 7.44 and for the geometric parameter Y=31.5 
chosen. Such a behaviour has also been observed by Mead[24] 
for the fundamental mode. Again, apart from the first re- 
sonant frequency, little change Is observed in the natural 
frequencies and loss factors between the fixed edge con- 
ditions and rbtationally restrained edge conditions. The 
effect of the boundary conditions is felt significantly only 
in the fundamental structural mode. 



Table 4.1 - Natural Frequencies and Loss Factors for Different Plate Boundary 

Conditions 
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Table 4.2 gives the variation of the natural frequen- 
cies and loss factors with the core shear parameter g* for 
simply supported edge conditions at x = + a/2. The other 
core parameters are kept constant. It is seen that it has 
the most significant effect on the natural frequencies and 
loss factors . The resonant frequencies increase with in- 
crease in g* in all the modes, the increase in frequency 
becoming less for high values of g* . It is so because, 
the core becomes stiff with increasing values of g * . The 
loss factor decreases with the modal number for g *= 1.0, 
while it increases with modal number for g* = 100.0 and 
1000.0. For g* = 7.443, the loss factor increases upto a 
certain modal number and then decreases. It is also obser- 
ved from the table that for a particular mode, there is an 
optimum value of g* for maximum loss factor. This optimum 
value increases for higher modes. Since the modal loss fac- 
tor depends on a complex fashion on the modal number and 
the shear parameter as can be seen from the expression for 
the simply supported edge conditions (equation (2.22)) no 
special reason can be attributed to such a behaviour. All 
that can be said is that damping of a sandwich plate is 
known to be mode dependent. Evidently when g* is changed 
the modes of the sandwich plate are changed and hence the 
variation of the loss factor with g * . 

Table 4.3 shows the variation of the natural frequencies 
and the associated loss factors of the sandwich panel in the 
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symmetric modes for the edges x = + a/2 fixed. The 
variation of the natural frequencies is very similar to 
that for the edges x = + a/2 simply supported. The loss 
factor is less in the fixed edges case than simply step- 
ported case for g*= 1.0, but they are greater for g * = 100. 
For the value of g* = 7.443, the loss factor for fixed end 
conditions is more than the corresponding value for simply 
supported conditions in the first fundamental mode, but 
for other modes the loss factors are less. Again it is very 
difficult to offer a physical argument for such a beha- 
viour, except to conclude that the modal loss factor is 
sensitive to the variations in g* and the form of the 
damped normal modes. 

Tables 4.4 and 4.5 show the variation of the resonant 
frequencies and loss factors with core loss factor, for the 
y wise edges simply supported and fixed respectively. In 
both the cases, increasing p has only a marginal effect 
on the resonant frequencies, but the modal loss factors 
increase almost linearly with pas is the case with the 
infinite sandwich panel. 

Tables 4.6 and 4.7 show the variation of the resonant 
frequencies of the sandwich panel and the associated modal 
loss factors with the core geometric parameter Y, for 



Table 4.4 - Natural Frequencies and Loss Factors for Different Core Loss Factors 
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Table 4,5 - Natural Frequencies and Loss Factors for Different Coro Loss Factors 

(x * + a/2, Fixed, g* = 7.443, Y = 31.5 ) 
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Table 4.7 - Natural Frequencies and Loss Factors for Different Geometric Parameters 

(x = i a/2, Fixed, g* = 7.443, p « 0.3 ) 
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the edges at x = + a/2 simply supported and fixed respecti- 
vely. It is seen that both the resonant frequencies and the 
nodal loss factors increase with increasing value of Y in 
both the cases . The nodal loss factors also increase with 
Y , which is expected because of the increased thickness of 
the damping layer. 

The cavity modal frequencies for the dimensions of the 
cavity are shown in Table 4.8. 

4.5.3 Comparison of Results in the Two Formulation 

The acoustoelastic formulation of the noise transmission 
and structural response problem represents a more exact ana- 
lysis of the coupled cavity-structural system. Hence, only 
a few results following the approximate foimulation given 
in section 4.2 are presented for purposes of comparison. 

Figure . 4 .3 shows the variation of noise reduction at 
the centre of the plate, with respect to frequency in both 
the Vaicaitis (approximate) and the acousto elastic formu- 
lations for the edges at x = + a/2 simply supported, and 
for an external white noise excitation. The same number 
of cavity modes and structural modes are included in both 
the analyses. It is observed that neglecting the effect of 
the cavity pressure on the panel motion does not alter 
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significantly the location of the dips in the noise reduc- 
tion curves which occur very nearly at the cavity res onant 
frequencies . At low frequencies, the approximate formula- 
tion overestimates the noise reduction values by 1 to 3 dB 
at the cavity resonances, while at higher frequencies it 
underestimates the noise reduction values by the same or- 
der. However , there is significant difference in the n on- 
resonant noise transmission characteristics of the plate in 
the two formulations. 

The structural response for the same conditions as 
above is shown in Figure 4.4. In the low frequency ranges 
(0 - 200 Hz) there is practically no difference in the 
structural response in the two formulations. In the higher 
frequency regions, the in-vacuo structural response obtained 
by neglecting the cavity pressure on the plate motion is 
more than the structural response obtained when the full 
acoustic-structural coupling is taken into account. This 
may be attributed to the fact that the cavity modes in cer- 
tain frequency bands act as dynamic absorbers for the struc- 
tural motion. 

The computational efforts involved in the coupled acous- 
toelastic formulation turned out to be even less than for 
the approximate formulation for the total number of modal 
summations considered. For the simply supported case the 
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CPU tine on a DBCIO computer for the acoustoelastic formu- 
lation considering a total number of 25 structural and ca- 
vity nodes turned out to be 30 sec . for 40 frequency cal- 
culations. In the approximate formulation it turned out 
to be of the order of 1 minute. 

Henceforth, the results obtained by the acoustoelastic 
formulation are presented. 

4.5.4 Effect of Plate Boundary Conditions 

Figure 4.5 shows the variation of the noise reduction 
at the centre of the plate for different panel boundary con- 
ditions at x = + a/2 and for a harmonic external pressure 
excitation with normal incidence. The trend of variation 
is the same for all the boundary conditions and beyond the 
fundamental structural resonant frequency, the noise reduc- 
tion values do not differ by more than 4 dB at any frequency 
for the different boundary conditions . Below the first ca- 
vity resonance, (nonzero frequency) dips in the noise reduc- 
tion curves occur at the first fundamental structural re- 
sonant frequency. The dip is most predominant in the case 
of transversely restrained edges and least noticeable for 
the case of fixed edges. The difference in the noise reduc- 
tion values corresponding to these two curves is as high as 
20 dB around these frequencies. The low noise reduction 
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value for the case of transversely restrained edges can be 
attribute! tc the low value of the modal loss factor as 
seen from Table 4.1. For the case of fixed edge conditions 
the high noise reduction value can likewise be related to 
the high modal loss factor corresponding to the first reso- 
nant frequency. For the case of rotat ionally restrained 
edges, there is virtually no change in the noise reduction 
curves between riveted and unriveted edges for the torsional 
stiffness parameter chosen. Hence they are indistingui- 
shable in the curves. The reason for such a behaviour can 
be found again in Table 4.1, where practically no difference 
is discernible in the natural frequencies* and loss factors • 
values for the two cases. 

The noise reduction curves for a random external pres- 
sure excitation and for different y wise edge conditions 
are shown in Figure 4.6. The nature of the curves is very 
similar to that obtained in the case of harmonic excitation, 
as the random pressure field is weighted uniformly in all 
the frequency ranges being a white noise excitation, and 
spatially uniform. The deflection power spectral densities 
at the centre of the panel for these boundary conditions and 
external random excitation are shown in Figure 4.7. For 
all the boundary conditions, the displacement spectral den- 
sity peaks nearly at the first structural resonant frequency, 
indicating that for the cavity depth adopted, the structural 






mply supported 
Nationally restrained 



OUNDARY CONDITIONS ON STRUCTURAL RESPONSE 



129 


resonances are not affected by the coupling action. The 
damping induced by the viscoelastic layer, flattens the 
peaks corresponding to other structural resonances. JL 
careful scrutiny of the results reveals that the displace- 
ment maxima occur very nearly at the in-vacuo structural 
resonant frequencies and if there have been shifts in these 
peaks due to coupling action, they are not noticeable. The 
displacement response for the transversely restrained con- 
ditions is generally more than for other boundary con- 
ditions upto a frequency range of about 600 Hz, persumably 
because of the low modal loss factors associated with this 
boundary condition. For high frequencies above 800 Hz the 
response spectral density becomes almost independent of the 
nature of boundary conditions. 

4.5.5 Effect of Core Loss Factor 

The change in noise transmission characteristics of 
the sandwich panel with respect to the core loss factor p 
is shown in Figures 4.8 and 4.9 , respectively for sim- 
ply supported and fixed edge conditions at the y wise edges 
and for external white noise excitation. For p = 0 re- 
presenting an elastic core, apart from the dips at the ca- 
vity resonances, dips corresponding to structural resonances 
also occur in the noise reduction curves. These dips are not 
present for p = 0.3 and 1.0, except for the first structural 



SMI 




132 


resonance. Even in that case, the noise reduction values 
are considerably higher for higher values of p than for 
p = 0. Increas ing p has little effect on noise reduc- 
tion values at cavity resonances. This observation agrees 
with that of Guy and Bhatt ach ary a [53], who in their ana- 
lysis of sound transmission through a cavity backed finite 
panel showed that the noise reduction at cavity resonances 
is independent of panel impedance. Figure 4.8 shows nega- 
tive noise reduction values for p = 0 around frequencies 
of 630 Hz and 650 Hz. It is so, because of the strong 
coupling of the (m = 5, n = l) structural mode which has a 
natural frequency of 628 Hz with the cavity mode (i=0, j=0,k=3 ) 
which has a frequency of 649 Hz. Negative noise reductions 
have also been reported by Guy and Bhatt acharya [53], and 
Mead [90] in connection with noise transmission studies of 
finite homogeneous panels. With increased p the noise re- 
duction is slightly improved. 

The displacement spectral density functions corres- 
ponding to the above two cases are shown in Figures 4.10 and 
4.11. Increase in core loss factor significantly reduces 
the structural resonant responses as can be seen from the 
figures. Corresponding to a frequency of 630 Hz a large 
structural response restats for the simply supported con- 
ditions because of the strong coupling between the struc- 
tural mode and cavity modes earlier referred to, for an 



Eq 0 1= 0 02 d 3 0 7 62 
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elastic core. With increasing values of § the response 
at this coupled frequency is considerably reduced as seen 
from Figure 4.10. 

4*5.6 gffect of Shear Parameter 

Figures 4.12 and 4.13 present noise reduction curves 
for different values of the core shear parameter g* , res- 
pectively for simply supported and fixed edge conditions 
at x * + a/2. The external excitation is a white noise 
excitation. For the values of g* chosen there are dips 
in the noise reduction curves corresponding to the first 
resonant frequency of the structure, in both cases of edge 
conditions. For g* = 100 and g* = 1000, there are other 
dips corresponding to structural resonances. This is due 
to the fact that the modal loss factors corresponding to 
those values of g* are very less. At low frequencies 
corresponding to g*= 1000 , negative noise reduction values 
are obtained around 126 Hz which is the first resonant fre- 
quency of the structure. It agrees with the earlier obser- 

*2 

vat ion in section 4.4 that for low values of and 
the noise reduction could have negative values as seen 
from equation (4.69). For g* = 1.0 the noise reduction 
at 225 Hz is very nearly zero indicating a strong empling 
between a structural mode ( ^3 = 228 Hz) and a cavity 
mode ( <j 3 qo1 = 216.5 Hz). For g* = 7.44 such coupling arises 
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between structural modes { _ 62? Hz, = 665 Hz ) 

and the cavity modes ( Wqo3 = 649 Hz and u> q2q = 660 Hz) 
which is manifested by nearly zero noise reduction cor- 
responding to a frequency of 650 Hz. In almost all cases 
considered there is a predominant dip around 650 Hz which 
is due to the coupling between the two proximate acoustic 
modes of OQ ^ = 649 Hz and w Q 2 0 = 660 Hz. The noise re- 

duction values are considerably higher for higher values 
of g* till about 125 Hz. In the fixed end conditions, 
there are coupled resonance conditions around frequencies 
of 200 — 225 Hz and 650 Hz as seen from Figure 4.13. 

Figures 4.14 and 4.15 respectively, show the response 
spectral densities for the above cases. The frequencies 
at which the response peaks are slightly shifted to the 
right from the in-vacuo resonant frequencies for the lower 
order modes. This is because of the proximity , of these 
modes. to the rigid cavity mode w 000 = 0. The closer the 
structural resonant frequencies are to & 000 * "the greater 
is the shift to the right. The structural response is 
generally lower for g* = 7.44 than for g*= 1 or 100 
because of the relatively higher modal loss factors, in 
the simply supported conditions. The structural response 
is the highest for g* = 1000. 






EFFECT OF SHEAR PARAMETER ON STRUCTURAL RESPONSE 
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4.5.7 Effect of Geometric Parameter 

Figures 4.16 and 4.17 show the effect of the core geo- 
metric parameter Y on the noise reduction inside the ca- 
vity. Three representative values of Y corresponding to 
thin, medium and thick cores are chosen. Figure 4.16 shows 
for the simply supported conditions, negative or zero noise 
reduction values for Y = 3.63, corresponding to frequencies 
—of 45 Hz, 225 Hz and 650 Hz. It is attributed again to the 
strong coupling between structural resonances and cavity 
resonances ( = 38.21 Hz, w OOQ = 0); (« * 1 = 206 Hz 

and co oq1 = 216 Hz) and (<*>*55 = 642 Hz and oi QQ ^ = 649 Hz, 

^2^ = 660 Hz) because of their close proximity. Such 

ozo 

coupling exists for the fixed edge conditions because of 

the close proximity of the structural modes ( = 624 Hz, 

* 

= 658 Hz) and the cavity resonances ( <a Q0 ^ = 649 Hz 
and = 666 Hz). In fact, in all the noise reduction 

ole 

curves, a clear dip Is shown at a frequency of 650 Hz 
because of the very strong coupling between the acoustic 
cavity modes = 649 Hz and ^qq - 666 Hz which are 

independent of the plate boundary conditions, damping and 
stiffness and depend only on the dimensions of the cavity. 
For very low frequencies below the fundamental resonance 
of the panel the noise reduction increases with increase 
in Y, but the effect of the ge cmetric parameter on the 
noise reduction at high frequencies is minimal. 
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Figures 4.18 and 4.19 present the structural response 
variation with Y for the two boundary conditions. For 
Y = 3.63 the first peak in the structural response is per- 
ceptibly shifted to the right by about 10 Hz both in the 
simply supported and fixed cases. It is due to the effect 
of the rigid cavity mode w 000 = 0 which adds to the stiff- 
ness of the plate. The coupled structural -cavity natural 
frequency can be approximately cal culat ed from equation 
(4.71) . The calculations agree well with the frequency 
of the first structural response peak. 

4.5.8 Effect of Cavity Damning 

Figures 4.20 to 4.23 show the effect of cavity damp- 
ing on noise reduction and structural response respectively. 
The other conditions are shown in the figures. Increase in 
cavity damping flattens the noise reduction curves at the 
cavity resonances and increases the noise reduction consi- 
derably at the lower order cavity resonances. For higher 
frequencies, the effect of cavity damping on noise reduc- 
tion is not appreciable. However, Increased cavity damping 
brings down the maximum noise reduction values as the cavity 
impedances at these frequencies increase with increased 
cavity damping. The peaks in the noise reduction curves are 
termed as ’cavity off resonances [53] in the literature 
and occur when the cavity modal Impedance attains a mi nim um 





FIG. 4.19 EFFECT OF GEOMETRIC PARAMETER ON STRUCTURAL RESPONSE 
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value. For higher frequencies the cavity damping does 
not have significant effect on noise reduction, because 
the modal damping coefficient defined in equation (4.11) 
varies inversely as the square of the cavity resonant fre- 
quency. Cavity damping has little effect on the structural 
response as depicted by Figure 4.21, which also agrees the 
result of Guy and Bhattacharya [533. 

4.5.9 Effect of Incidence Angle 

The noise reduction curves shown for harmonic excita- 
tion so far correspond to normal incidence of the external 
pressure excitation. Figure 4.24 shows the noise reduction 
curves for different incidence angles. It is observed that 
the incident angle has practically no effect on the noise 
reduction upto a frequency of 650 Hz. Above that frequency 
the noise reduction values increase with increasing angles 
of incidence which behaviour is quite different from that 
of the infinite sandwich panel. However the dips in the 
noise reduction curves which may be likened to the coin- 
cidence of finite panels described by Bhattacharya, Guy 
and Crocker [54] occur at frequencies independent of the 
incident angle. 




152 



000 


153 


4.5.10 Effect of Cavity Depth 

The effects of changing the cavity depth on the noise 
reduction and structural response are shown in Figures 4.25 
and 4. 26 respectively. Increasing the cavity depth adds to 
the number of cavity modes and decreasing the cavity depth 
decreases the number of cavity modes in the frequency range 
of Interest . Accordingly the dips in the noise reduction 
curves change depending on the coupling between a structural 
mode and cavity mode. The change in the cavity depth has 
negligible effect on the st ructural response for the depths 
shown. However, for d = 1.524 m the coupled natural fre- 
quency of the cavity panel system is shifted to 75 Hz, because 
the proximate acoustic cavity mode of 108 Hz is greater than 
the fundamental structural frequency of 81 Hz and hence adds 
to the total inertia of the panel . ,| 

4.6 CONCLUSIONS 

The noise transmission characteristics and the struc— : 
tural response of a finite sandwich panel with constrained 
damping layer, backed by a rectangular cavity is invest!- j 

gated in this chapter in an acoust oelastic formulation | 

which takes into account the coupling between the inter- | 

nal sound pressure field inside the cavity and the struc- ! 

tural motion. An exhaustive parametric study with respect j 
to the core parameters is made. j 
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The noise transmission characteristics of finite 
sandwich panels backed by a cavity are very similar to 
those of finite homogeneous panels. But because of the 
presence of the damping layer the noise transmission cha- 
racteristics of the sandwich panel are considerably im- 
proved at the lower order structural resonant frequencies, 
especially at the panel’s fundamental resonance. A pro- 
per choice of the core shear parameter g* ensures large 
modal loss factors which improve the noise transmission 
characteristics of the sandwich and reduce the structural 
response at structural resonant frequencies. Thus the 
analysis presented indicates the existence of an optimum 
value of B* for maximum damping effectiveness correspond- 
ing to the different structural modes. Increase in the 
core loss factor , has only marginal effect on the sound 
transmission at higher frequencies. But the effect of f. 
is strongly exhibited in the structural response, which 
is considerably reduced by increase in P at the struo- 

tnraX resonant frequencies • 


The effect of the plate boundary conditions is no- 
ticeable only in the lower frequency ranges and is mar- 
ginal at higher frequencies. This chapter clearly shows 
the efficacy of using the damped normal modes in response 
calculations and noise reduction computations, to take 
l of the arbitrary boundary conditions and arbitrary 


damping* 
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This chapter also presents a simple matrix inver- 
sion scheme to solve the coupled acoustoelastic problem, 
both for harmonic and random excitation which avoids the 
necessity of carrying out an eigenvalue analysis. 
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shells by including the effects of external flow, reflec- 
tion and scattering of incident wave, internal cabin pres- 
surization, cavity resonances and structural damping. 

A survey of the literature, reveals that there is no 
work on the sound transmission characteristics of sandwich 
shells with additive damping treatments. In continuation 
of the investigation on the scund transmission throu^i in- 
finite sandwich panels am finite sandwich panels conside- 
red in Chapter 3 and Chapter 4, the noise transmission and 
structural response of finite closed circular cylindrical 
sandwich shells are analysed in this chapter. In addition 
to the treatment of the sandwich shell with constrained vis- 
coelastic damping layer, the case of a layered shell with 
unconstrained damping treatments in different configurations 
of the damping layers is also considered in this chapter. 

v ui i interaction between the internal cavity pressure 
field and the structural motion has been taken care of in 
an acoustoclaat ic formulation. Only the axisymmetric vibra- 
tions of the shell are considered in the analysis as the 
variation of the external pressure field along the circum- 
ference is assumed to be -uniform. The shell ends are 
assumed to be simply supported which facilitates closed 
form solutions of the natural frequencies and the modal loss 
factors of the layered shell. There is no restriction on 
the formulation, if other boundary conditions are specified 
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at the shell ends, except that it would involve iterative 
procedures similar to that outlined in Chapter 2, in the 
determination of the natural frequencies and damped normal 
modes. Results are presented as noise reduction curves 
for both harmonic and random excitation. The shell response 
is also plotted as a function of frequency. The frequency 
range of interest has been restricted to 1000 Hz. 

5.2 EQUATIONS OH MOTION FOR LAYERED SH5LL WITH UNCONSTRAINED 
DAMPING LAYER. TREATMENT 

The layered cylindrical shell considered in this chap- 
ter consists of a thin circular elastic shell of finite 
length 1, to which is applied an unconstrained viscoelastic 
damping layer. Three different types of layer construc- 
tions arc- considered. 1. Outside coating (OC), the damping 
treatment being applied on the outside of the elastic shell, 
(Fig. 5. 1(a)}. 2. Inside coating (IC), the damping treatment j 
being applied on the inside of the elastic shell, (Figure 5»l( 
3. Two-side coating (TC), the damping treatment being sym- ■ 
metrically applied both inside and outside of the elastic 
shell, (Figure 5.1(c)). j 

The equations of motion of the layered shell in axi- ; 
symmetric vibration have been adopted from Markus [37], I 

who derived them under the following assumptions. j 



Median srface 



FIG. 5.1 DIMENSIONS AND COORDINATE SYSTEM OF LAYERED SHELL a) OUTSIDE COATING (OC 
t>) INSIDE COATING (1C) c) TWO SIDE COATING (TO 
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(i) The thickness of the shell h^ is much smaller than 

•pi 

the radius of the median surface r, that is, (h^/r) « 1. 

(ii) The viscoelastic material has a complex modulus 

** 3* (1+ip ), where p is the loss factor of the visco- 
elastic material and S 2 is the storage Young’s modulus 
of the damping material. 

($2 and p are assumed to be frequency independent in the 

analysis ) . 

(iii) The face layer to which the damping treatment is 
applied is assumed to be elastic and non-d iss ipat ive . 

(iv) The cylinder motion in the transverse direction for 
all the layers is the same. 

Considering the displacement continuity at the inter- 
faces of the layers and the equilibrium of forces and mo- 
ments, and eliminating the longitudinal displacement the 
following equations for the shell motion in the transverse 
direction for the three cases of dampixg layer configuration 
OC, IC and TC can be obtained. 

A w 17 - (B/r) w 11 + (C/r 2 ) w + Jiw = p(x,t) (5.1) (00 ) 

A w 17 + (B/r) w IX + (C/r 2 ) w + pw = p(x,t) (5.2) (IC) 

D w 17 + (C/r 2 )w IX + pw = p(x,t) (5.3) 

where, A,B,C and D are constants depending on the extensions] 
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rigidities of the different layers in the composite shell, 
their Poissons ratios and thicknesses, p is the surface 
mass density of the composite shell. 

The constants A, B, C and D arc given by, 

1 O ‘ 1 o o 

A = j Ojh^ + gr oo, h! + d (oc^ + cc,) - d( ctjh^ - ) 

(5.5) 

B SB £ a x( V 2“ V l)( a 2* b i +a: 2 h 2)/ ( a 1 +a 2*) 

+ 1 (a^Vg - a 1 b 1 v 1 ) + d(a.jV^ + a^) 


+• do^ ( v 2 - v 1 )( a 2 - a 2 V (“]_ + a 2^ (5.6) 

C = a 1 + a 2 - (a^v-j + a 2 v 2^ a l V l 4 ‘ a 2 V 2^/^ a l +0 ^ ) (5-7) 

D = 12 (“l 31 ! + “2^2 + 3a 2 h i ^ + (5.8) 

where, d - \ (a^ - a* h 2 )/(a 1 +- ) 

a I = S i h i /(1 ” v i ) ( 5-9) 

a 2 = E* h 2 / (1- n§) (5.10) 


and E-^ Is the Young’s modulus , v^, the Poisson’s ratio 

JL"U 

and h. the thickness of the i ^ layer. 

I 

For an external loading of the form p{x,t )=-iT}pw f it 
can be shown that the motion of the plate is harmonic and 
can be expressed in separable form 

i<i> t 


w = ¥(x) e 


(5.11) 
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Substituting equation (5.11) into equations (5.1) to 
(5.3), one obtains , 

k - (B/r) W 11 + (C/r 2 )¥-|A w* 2 (l+iq) ¥ = 0 (5.12) 

k Y 17 + (B/r) W 11 +• (C/r 2 )W-ji M * 2 (l+iq) V = 0 (5.13) 

D ;: IV + (C/r 2 )y - pui* 2 (1+iTj) V = 0 (5.14) 

where w* is the square of the natural frequency and 
q the associate loss fact or. 

The natural frequencies , loss factors and the damped 
normal modes would depend on the boundary condit ions at 
the ends of the shell . As the main interest pertains to 
the noise transmission through the shell into the cavity, 
and as it had been observed in the case of finite plates, 
that the behaviour of the noise transmission does not vary 
qualitatively with the plate edge conditions , the analysis 
is restricted to simply supported end conditions . For such 
a case, the natural frequencies and the loss factors can 
be obtained by closed form expressions. It can be shown 
that for simply supported ends the damped normal modes take 
the form 


^.(x) = ^ sin 


E-ttX 


(5.15) 


for all the three cases of layer constructions OC,IC and TC 



165 


Substituting equation (5.15) into equations (5.12) to (5 .14) 
the natural frequencies and the nodal loss factors are given 
by 

^ <!«*) - CA-y + (B/,) (5-16) (OC, 

u* 2 (1+1%) = [A^f- +{B/r) S-f- +(C/r 2 )]/|i (5.17) (IC) 

1 1 

% 2 (l+irfe) = [D +(C/r 2 )] (5.18) (TO) 

The natural frequencies and associated loss factors 
for various elastic and geometric properties of the layer 
materials are given in Tables 5.1 to 5.3 . 

5-3 EQUATIONS OF MOTION 0? DAMPED SANDWICH SHELL 

The sandwich shell considered is a thin circular cylin- 
drical shell (Figure 5 . 2 ), consisting of three layers , two 
symmetric face layers of equal thicknesses and same material, 
with a constrained viscoelastic damping layer. The equa- 
tions of motion of the sandwich shell in axisynmetric vibra- 
tions have been derived by Pan [36] following an analysis 
of DiTaranto [12] for three layered sandwich beams. The I 

following assumptions were made in the derivation. 

i) The shell is thin and hence the radial displacements 




taDle 5.1 - Resonant Frequencies and loss Factors of layered Shell 

(1^/h, = 0.1, p = 0.3, e£ = 4.14 x 10 6 ti/m 2 ) 
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Table 5.3 - Resonant Frequencies and loss Factors of layered Shell. 0.0. 
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for all the three layers are assumed to he the sane and 

6QU&1 W* 

ii) The bending stresses in the face layer are much 
larger compared to the direct bending stresses in the 
viscoelastic core. 

iii) Displacement continuity is maintained at the inter- 

« 

faces of the three layers . 

iv) Inertia forces in the longitudinal direction are 
negligible* 

Neglecting the ertcnsional force carried by the vis 
coelaatic layer and neglecting the longitudinal inertia, 
force equilibrium in the axial direction results in 


("i 


B^h^ 


v i!> + ?^S- V 3l ) = 0 


( 5 . 18 ) 


th 


E l h l 

where, E^lu and u ± refer respectively to the Toung 
modulus ," 1 "thi ckness and mid-surface displaoaaent cf the i 
layer respectively. A prime indicates diff erentiationwith 
respect to r and % is the Poisson's ratio of the i 
layer, w is the radial displacement of the shell and 
a is the radius of the midsurface of the core layer. 


The displacement 
sed by the relation 


continuity at the interfaces is expres- 
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<9 + *) 


( 5 . 19 ) 


where , <p = Sw/ ?x is the rotation of the section, assumed 
same for the layer 1 and 3 and ? is the shear strain of 
the viscoelastic core layer. 


Considering the equilibrium of forces in the longitu- 
dinal direct ion cf an element of layer 1 it can be shown 
that (Figure 5.3) 

* v “11 ) ( 5 . 20 ) 

Q j.— vj ^ B. 

where G = G* (1 + ip ) is the complex shear modulus of the 
viscoelastic core . 

Making use of the fact that 3^ = = B and v 1 = = v 

and using the continuity condition (5 .19 )in conjunction 
with equations (5.18) and (5.20) one gets. 


* bA 

6(1- v 2 ) 


I t X 


>u£ 


+ 2 


v w 
a 


Bh-ho 

+c(hl+h2) " ^ b 7) lw ' ,=0 


(5.21) 


The bending moment M x at any section is obtained by 
taking moments of the axial forces and is given by, ( assum- 
ing extensions! forces in the core damping layer to be zero) 


M x = 


B iu 

1(^+1^ )(u£ - — )-D t W 


i t 


1-V-. 


a 


( 5 . 22 ) 
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where , 


D, = 


E h£ 


6(1- v ) 


(5.25) 


For the axisymnetric vibration of the shell, the equa- 


tion of motion can be represented as 

K i' + I + 


(5.24) 


W ir “hr> 1 


™ferent.i 


of the whole sandwich given by 


■>rce per unit length 


M = - 2 E ( w/a) 


(5.25) 


p. is the surface mass density of the composite shell and 
p(x,a,t ) is the external pressure loading on the surface 
of the shell. 


Substituting equations (5.22) and (5.25) in equation 
(5.24) and using equation (5.21) to eliminate from the 
equation, the equation of motion in w can be written as 
2G(l- v 2 )w iV Eh- (h.+h-) 2 

W VX [1+ * ] 

E h^ 2B t (l- v^) 

B h. tt G(l- v 2 ) b- 

+• — £ W - —X ^ w 


2a D^. a h^hg 

•M. «V 4 2 w, l 

m m &*mm**mm**§' - nr-T.j-- j -r i i n n. ***** mmmmm — » mmm. mmmmmi | 






1 -2 

i_ Ue. 

B. 0x 2 


(5.26) 
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Equation (5*26) can be rewritten as 


- g‘ ( 1+Y )W 1V + 


i_^) ,, isg’d-v 2 ) 


w - 


** e 2 w' * 
t o~ 

D t at* 


.2 1.2 

i»i -M- y 1 4» ~ .. .ip. 

at 2 u e x 2 


(5.27) 


where g 1 


^1^2 


is the shear parameter and 


2 

Y = 3 (l+ c— ) is the geometric paraneter which 

n l 

are the sane expressions as for the sandwich plate. 

For an external loading of the form p(x,t)=in|i 

2t 

it can be shown that the plate motion is harmonic and can 
be expressed in separable form in damped normal modes as. 


w = ¥(x) e 


i <*>t 


(5.28) 


Substituting equation (5.28) into equation (5.27) one 


obtains, 


V 71 - g* (1+Y )W^ 7 — ( V^-g’ W)[— — (1+iri ) - 


12(1- v^) 


■] = 0 


(5.29) 

where w* is the square of the natural frequency and q 
the associated loss factor. Equation (5.29) is of similar 

of motion for sandwich beams with con- 
i’ [13]. 
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The natural frequencies, the loss factors for simply 
supported end conditions can he obtained in closed fora 
by substituting pure sinusoidal modes of the form 

mux 

¥ m (x) = Swj sin — (5.30) 

into equation (5.29). The resulting expression for the 

"til 

complex natural frequency in the m mode is given by 

(fS) 6 + g.(l+Y)(fE) 4 + l|U- v .. !) { ( aE)2 + g.] 

a h^ 

( f *) 2 + g . f 

The natural frequencies and the associated loss factors of 
the sandwich shell for typical values of p, g and Y are 
presented in Tables 5.4 to 5.6. 

5.4 ACOUSTOELASTIC FORMULATION 
5.4.1 Acoustic Solution 

Consider the closed cylindrical cavity shown in 

2 

Figure 5.4, occupying the cylindrical volume ¥ = na 1. 

It is assumed that the wall formed by the layered shell 
or the sandwich shell is flexible and the two ends of 
the cylindrical are closed with rigid walls. The pres- 
sure p c inside the cylindrical cavity is governed by the 
acoustic wave equation in cylindrical polar coordinates. 


2 (1« v = ^ 



Table 5,4 - Resonant Frequencies and Loss Factors for Different Shear Parameters 

(Sandwich Shell ) 


175 


i|i > 


I 


1 

1 











1 


1 


a> 


kx 

*y 

CM 

H 

H 

H 

ri 

H 





I 


i 

1 

i 

1 

i 

1 

I 

I 

1 


s 


I 


o 

O 

o 

o 

O 

O 

O 

O 

O' 

O 





1 

1 

** 

H 

X 

a 

d 

H 

K 

A 


H 

« 

A 

H 

K 

A 

© 


1 



1 


H 

K\ 

£> 

OX 



P 

KX 

OX 

I 



1 


kx 

kx 

in 

cc 

co 

CM 

QX 

MD 

KX 

O 




! 


KX 

CM 


CM 

OX 

CM 

KX 

in 

c*- 

OX 

H 

I 



I 

I 

I 

! 


• 

• 

* 

• 

• 

♦ 

» 

• 

• 

• 


1 


o 













1 


o 

I 














I 

o 

1 

m 










O 


! 

in 

! 

Dd 

H 

o 

CM 

CM 

H 

Ox 

ox 

n- 


O 




1 


H 

H 

X0 

CM 

KX 

in 

H 


H 

in 

in 

§ 



1 


* 

m 

• 

* 

• 

• 

• 

• 

« 

• 

m 




I 

a 

CM 

t <\ 

ox 

OX 

OX' 

KX, 

H 

CO 

H 

in 

€Q 

1 



1 

I 


O 

o 

o 

CM 

MO 

KX 

CM 

CM 

in 

CD 

CM 






<- 


■sK 


in 

VO 

f> 

ID 

ox 

H 

1 

1 


1 

I 

I 












H 

1 



I 


in 

KX 

CM 

CM 

H 

H 

H 

H 

H 

H 

H 

1 




! 1 

I 

1 

I 

I 

I 

i 

1 

i 

1 

1 

1 

in 


1 


O 

O 

o 

O 

o 

O 

o 

O 

O 

O 

O 


* i 

H 


S 

I 

a 

i H 
X 

A 

$ 

3 


d 

$ 

A 

3 

d 

K 

r 

K X f 


? 

it 

[ o 

CXI 

VO 

P 

r- 

! Ox 

s> 

n- 

r-4 

in 

in 

i 



« 


CXI 

in 

H- 

CM 

OX 

o 

in 

KX 

0X 





1 

I 


i H 

kx 

C\1 

'*>0 

H 

H 

CM 

KX 

^K 


in 

l 

11 

m 


* 

* 

m 

* 

# 

# 

• 

* 

# 

w 

* 


H 

I 


! 











! 


Hi 

1 
















I 


i 











I 



1 

& 

i A 

71 

88 

ox 

CM 

CO 

H 

in 

in 

CD 

in 




f 

w 

CD 

O 

MO 


!> 

CM 

KX 

in 

I 



I 


• 

• 

• 

• 

• 

* 

• 

• 

* 

m 

* 


k\ s 


1 

„ a 

f CM 

CM 


CD 

in 

KX 

in 

O 

0 


KX 

I 

* 


I 


o 

o 

o 

O 

H 

CXI 

KX 

in 

r- 

OX 

CM 


o 1 

li i 

ax i 


1 

1 

t 

i 










in 

1 

i 


I 

* 

I 











1 

1 



1 


in 

kx 

to 

KX 

CM 

CM 

CM 

CM 

CM 

CM’ 

CM 



1 


s 


1 

i 

1 

I 

1 

I 

i 

i 

I 

1 

I 

I 



1 



0 

O 

O 

0 

O 

0 

O 

O 

0 

O 

O 

i 


1 

1 


f 

A 

•A 

tf 



A 

X 

A 

A 

A 

A 

1 


1 

1 

CO 

in 

CD 

© 

H- 

KX. 

H 

KX 

CM 

KX 

O 



1 




KX 

CM 

VD 


H 

Ox 

O 

•in 

CM 

© 

« 

0 

1 


1 


H 


£0 

H 

CM 

CM 

KX 


tn 

in 

H 

I 



• 

* 

• 

m 

* 

* 

♦ 

* 

• 

* 

* 

• 


s 

1 


I 












1 


1 

****** 

f 














1 

m 

0 

KX 


KX 



KX 

VO 

0 

H 

CM 

1 


1 

« 

1 

H 

CM 

MD 

XO 

VO 

Cxi 

H 

0 

co 

H 

VO 



i 



• 

* 

♦ 

• 

• 

* 

* 

* 

w 

• 

* 

I 


1 

1 

a 

I 

CM 

CM 

CM 

KX, 

'in 

ox 

in 


XO 

H” 

\o 





O 

O 

O 

O 

O 

O 

H 

CM 

KX 

in 

m*' 

l 

I 

l 


1 

i 


« 

*sf* 



"'st* 



xK 








1 












I 

#bO 



1 












1 




I 












1 




1 












I 



a 

H 

KX 

in 

IH 

ox 

H 

KX 

in 

l> 

OX 

H 



1 






H 

H 

H 

H 

H 

cxi 



176 


ca 

U 

o 

+> 

o 

as 

Pi 

Q 

CO 

O 

H 

o 

U 

o 

o 

rt 

o 

Pi 

© 

<H 

•r-f 

Q 

Phi 

O 

© 

Pi 

O 

-P 

O 


I I 

i«ft 


I I 
O j . 
♦ 1 * 

H 1 1 

11 1 , 

CQ_ 1 « 

i a i 

i w 

i i 

! a, 

I *N * 
I 


A 

1 

O 

f~i 

K 

A 

A 

CM 


*? 

O 


A 

O 

A 


CM 

H 

H 

H 

1 

I 

1 

1 

o 

H 

O 

H 

O 

H 

O 

H 

H 

M3 

K 

A 

K 

ID 

K 

M3 


H 

H 

in. 

t- 

CM 

H 

M3 


T 1 

o 

K 

c~- 

H 

A 


CD 

'CXI 


H 


vo 

H 

H 

H 


O 

CD 

HI 


H 

H 

O 

CO 

CXJ 

A 

A 

X0 

o 

A 

03 

Oi 

M3 

OX 

'OX 

• 

CD 

<* 

m 

CD 

A 

m 

o 

a- 

• 

* 

CM 

o 

• 

CXI 

o 

-e- 

* 

in 

o 

■H 

* 

A 

o 

XD 

H 

<* 

CM 

*<c 

M3 

a 

H 

A 

H 

t- 

m 

0X 

e\i 

a 


1 


H 
H 

cj © 


I 

I 

! 

I 

^ I 

a i 

• 

H I 
A 


II 


A 

o 

II 

Wl. 




LT\ 

i 

o 


o 

CM 

H 


? 


cu 

I 

o 


CM 

I 


T 4, 

o 


T 1 

o 


CXJ 

A 


C- 

Ox 

A 

*#• 

CM 

A 

A 

CM 

* 

A 

* 

H 

• 

H 

«* 


-1 

H 

H 

H 

r 

i 

I 

i 

o 

o 

o 

o 

H 

X 

£“**■* 

3 

c- 

H 

H 

3 

A 

C— 

A 

A 

A 

CXJ 

A 

<- 



t 


o 

ITS 

"Ml* 

IA 

* 


I 

1 

I 

I 

I 

I 

I 


m 

m 

o 

H 


30 


I 


si 
o 
•H 
£ 

d £ 

d <3 

CO 


*rl 

O 

£ 

© 

3* 

C* 

© 

Pi 

pH 


a 

0 

m 

© 

Pi 

1 

a 

a 

© 

H 

■3 

£4 


I 

O 1 

* 

ox i 

H 

H I 


11 

♦3fi 


ts! 

W l 

H 

H 

F - 

88 

A 

CD 

CM 

O 

,68 

H 

• 

A 

* 

A 

CM 

• 

03 

A 

* 

A 

A 

• 

« a ! 

• 

CXJ 

O 

* 

CM 

O 

404. 

€0 

o 

A 

H 

<r 

A 

CM 

M* 

A 

A 

O 

A 

o 

t> 

A 

A 

CXI 

A 


I O 1 
J 


I 


O 1 P" I 

11 I 1 

CEL I I 

I 

I ^ I 
I N 

1 H ? 

i i 
i g 
1 *H 1 
1 

1 I 


O 

o 

o 

o 

o 

O 

O 

o 

♦ 

o 

# 

O 

• 

o 

* 

• 

o 

* 

o 

* 

o 

• 

o 

o 

o 

o 

o 

o 

o 

o 

H 

H 

o 

C“- 

A 

CO 

CO 

H 

A 

t- 

A 

A 

CM 

m 

A 

MD 

* 

CXJ 

H 

• 

.26 

m 

• 

CXJ 

o 

m 

CM 

o 

• 

o 

CD 

o 

<- 

H 

A 

CM 

A 

A 

O 

A 

kI* 

o 

H* 

& 

A 

CM 

A 


I A IA 


ox 


HI 

H 


A 

H 


H 


C- 

H 


ox 

H 


H 

CM 


I 

1 

I 

I 

1 

1 

1 

1 

1 

i 



Table 5.6 - Resonant Frequencies and Loss Factors for Different Geometric Parameters 

( Sandwich Shell ) 
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*2 

is. + 


(5.32) 


with the boundary conditions 


3 w(x, 8) 


at r = a 


(5.33) 


£ =0 at x = 0,1 

dx 

Equation (5-32) for homogeneous boundary conditions has 
normal mode solutions P^ e x ^ ' which satisfy the equations 


2 T? 

V F n 


- (to ,/c f F 
nA / n 


(5.34) 


where v is the Laplacian operator in cylindrical coor- 

jl-U 

dinates, is the n ^ cavity resonant frequency and 

the corresponding cavity mode. 


The cavity normal modes satisfy the orthogonality re- 
lation of the form 


F_ F dV = M , 6 


r n 


hA rn 


(5.35) 


where 6 rn represents the Krone cker delta. 


The acoustic pressure inside the cavity can be ex- 
panded similar to the case of the plates in terms of the 


acoustic normal modes F_ as 

XX 


PC £ 


P P 
n n 


(5.36) 
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The pressure coefficients P„ in the acoustic nor- 
* n 

nal node expansion are governed by the set of ordinary 
differential equations 


P +• 2t « , P + 
n ihx tux n 


2 « 
10. , P„ 
nA n 



(5.37) 


where. 


n 


Pc V 


/ 

¥ 


p„ d¥ 
n 


(5.38) 


W n = — / w F n dAj, (5.39) 

Ap Ap 


where Ap is the surface area of the flexible cylindrical 
wall and , is the nodal damping coefficient. 


5.4.2 Structural Solution 

The response of the cylindrical shell (layered and 
sandwich) can be expanded in the damped normal modes as 


w( x,a,t ) = 2 q(t ) <p. 

m 


n 


(5.40) 


where. 


<?■ 


m 


sin 


EHEX 


The normal coordinates are governed by the set of ordinary 
differential equations 




(5.41) 
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where, 

M E S - { 1*4 dA 

I’m dA 


Sae ~ “ *f 


i»y 


P e Tn 4A 


and p is the external pressure loading on the cylindrical 
shell. It has also been assumed in the analysis that the 
radiation loading on the cylindrical surface external to 
the cavity is negligible compared to the external pressure 
field. Substituting equation (5.40) In equations (5. 38) 
and (5.39) and equation (5.38) in equation (5.41) the 
coupled acoustoelastic equation in the pressure coeff ici ents 
and the damped normal coordinates are obtained as 


p + 2e 
n x 


* P 

. “nA P n + “nA P n = - f 

2 L €j 

ns in 
a 

(5.42) 

p °% 

PL ™ 

n nm qE 

M nA ^ 

(5.43) 


where. 


mu 


1 P n dA 


(5.44) 


"S' % 

are the acoust o-st ructural modal coupling coefficients 
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5.4.3 Modal Backing Coefficients ani Other Modal Parameters 

It can be shown that the acoustic normal node 7 , fcr 

n 

the cylindrical cavity ± 3 of the form 

? n (x,r, 0 )= cos cos j 0J^( Xj k r ) (5.45) 


where as seen from the equation, the index n includes in 

itself the three indices i,j ,k corresponding to the axial, 

circumferential and radial directions. In the foregoing 

equation J. represents the Bessel function of first kind 

"fell 

of order j and = am^/a where a^ k is the k root of 

the equation dJ\/dr = 0. 


Since only the a xi-sy metric vibrations of the shell 
is considered it is sufficient to admit in the acoustic 
mode only the j = 0 term in equation ( 5-45 ) as the contri- 
butions from the modes corresponding to j = 1,2, etc. to the 
generalized force or the coupling coefficient will vanish. 
Then the cavity modes can be expressed as 


F n (x,r) 


cos 


ittx 

1 


J ( 


X ok r) 


(5.46) 


The generalized mass of the structure in all the 
structural modes is given by 

2 7 | 1 

= / / \iam- Cin 2 ^ dx = ^aL (5.47) 



The cavity generalized mass from. equation (5.35) is, 


i a 2u 1 ~ p . 

: ^ = -5T l l ( TJ o ( \k r > cos ¥ SxJrie <5.48) 

%a 1 o o o 


which on integration yields, 

Jot *ok 3 > + 4 < *ok a > 


M nA - 


for i ^ 0 ani k ^ 0 


J o ( X ok a5+ ^ X o ^a) for 1 = 0 and k ^ 0 


for i ^ 0 and k = 0 
4 


(5. 


= -y for i = 0 and k = 0 


The cavity frequencies are given by 


a) 2 _ to 2 - c 2 F x 2 + (— } 2 ] 
nA “ ijk “ c 1 A jk + v 1 ; J 


(5.50) 


Ass xrain g the external pressure loading to be uniform 

i tot 

and harmonic of the form p e (x,a*t) = p Q e , the gene- 


ralized force 


2 m 1 2 si iwt 

Q jaB = / / P e (x,a,t ^(xjdxade = P 0 (l-cosmit )— e 


o o 


(5.51) 


The acoustic structural modal coupling coefficients 
L are given by 

TUB 


J o ' * 0 k a) x - 2 -* 


f f sin S3S cos ixade 


lux 


nm 


o o 


(5.52) 
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L 


m 


V X ok a) j 1 

2s + i 


COS (Erl-i) % 


+ 

m + i 

cos (n-i) s 
n - i 



( 5 . 53 ) 


The above equation has nonzero values only for the con- 
binations (i = even, n = odd) or (i = old, n = even) for 
whi ch ,equat ion ( 5 . 53 ) reduce s to 


L 


nn 


2a J o ( X ok a) 

% (n 2 - i 2 ) 


(5.54) 


Frees equation (5.51) it is seen that the generalized force 
Q w?i has contributions only for the odd order structural 
modes. In the solution of the acoustoelastic problem only 
the combinations (i = even, m = odd) cavity and structural 
modes are considered in the analysis . 


5.4.4 Acoustoelastic Problen Solution 

Equation (5.40) and (5.43) can be expressed in matrix 

# 

notation similar to that in Chapter 4 (equation 4 . 45 ). 


Hence, following the steps adopted in section (4.3.3) 
for a harmonic external pressure loading the cavity pressure 
can be determined as 


*p "P 

2 iw t v n x n 
pc e ^ 

n M nA 


(5.55) 
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and the structural response as 

v (x,a,8,t) = e iut | \ <(J*) (5.56) 

and the noise reduction as, 

HR = 20 log (5.57) 

where P^ and are the con pi ex amplitudes of the harmonic 
pressure coefficient P w and the structural generalized co- 
ordinate respectively. 


Similarly for a stationary uni font random external 
pressures with constant power spectral density function S , 
the power spectral density of the cavity pressure is given 

by 

I 

S (x,r,e,u>)= p 2 c 4 {a} [S pp (w)] [a*j (5.58) 

°^ c (lxn) (nxn) (nxl) 

The structural response power spectral density of the cylin- 
drical shell is given by 


s w (x #a ,e,») = {Y} [S qq ( «)] 

(1+m) (nxn) (nxl) 

where jaj. and | are vectors of orders n and m 
respectively with elements , 


(5.59) 


F n (x,r, 0) 


**ni 


and Y m = <p n (x f e,9) 
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and the matrices [S ?IJ (y)] and [S c ^( w )] are cross spec- 
tral density matrices of the pressure coefficients and 
structural coordinates, whose elements can be obtained 
from the cross -spectral density matrix of the generali- 
zed forces [S^qCoj )] in a similar fashion as explained in 
.. (mxn) 
section {4.3 .3). 

The corresponding cross acceptance function m ) 

between the ei n and r a structural node is given by 

u> )= / / r(x 1 , 8^ xg, e 2 i w 0 1 )q3 r (x 2 , ^Jd/^d/xg 

At? Ap 

* ( 5 . 60 ) 

where for r (x^, e^* Xg, 0 2 ; «) = 1.0 

, 2.2 

4 a 1 

I ( cu ) = (l-cosmu ) ( 1-cos TTt) (5.61) 

mr 

mn 

which has nonzero values only for m = odd and r = odd. 

The solution of the acoustoelastic problem can be 
effected by a matrix inversion procedure explained in 
Chapter 4. 
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5.5 results all discussion 

The material am geometric properties of the shell, 
both layered and sandwich, adopted in the analysis are 
given below. 

a) Layered Shell (Unconstrained damping treatment) 

Young's modulus of face layer B = 7.24x10^^1*/”^ 

Density of face layer material P = 2770 Kgm/m^ 

P 

Thickness of face layer h^ = C.002 m 

Mean radius of the shell r = a = 1.83 m 

Density of damping layer C{J = 0.2 

Loss factor of the damping layer 

material p = 0.3 

Length of the cylinder 1 = 2.0 m 


The other parameters like the thicknesses of the damping 
layers in the three configurations 0C,IC and TC, the 
storage Young's modulus of the damping layer are varied 
in the analysis. 

b) Sandwich Shell (Constrained damping treatment ) 

Young's modulus of face layer E = 7.24xl0^N/m^ 

*5 

Density of face layer material p = 2770 Kgm/m 

P 

Thickness of each face layer h^ = .001 m 

Mean radius of the shell a = 1.83 m 

Density of constrained damping layerp c = 0,2 p 


( core) 

Loss factor of the damping layer 

material 

Length of the cylinder 


p = 0.3 

1 = 2.0 m 
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The shear parameter g and the geometric parameters Y 
are varied in the analysis. The core loss factor p 
is also varied in the analysis. 

In both the cases, the cavity modal damping coefficient 
in most of the numerical investigations is taken to be 
e ool = The effect of cavity damping on the noise reduc- 

tion inside the cylinder is also studied by varying the ca- 
vity modal damping coefficient. 

5-5.1 Natural Frequencies and Loss Factors 
a) Layered Shell : Table 5.1 shows the natural frequencies 
and associated loss factors of the layered shell in the 
three configurations of 00, IC and TC of the unconstrained 
damping layer treatment. In the two sided coating the un- 
constrained damping layer is assumed to be symmetrically 
applied to the base layer on either side with damping layer 
thickness equal to half the thickness of the damping layer 
in the case of outside and inside coatings . The thickness 
of the damping layer is one tenth of the thickness of the 
face layer. The storage Young’s modulus of the damping 
layer is taken as = 4.14 x 10 N/m . 

It is seen from the table that the damping layer 
configuration does not have any effect on the natural fre- 
quencies and loss factors as they remain practically the 
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sane for all the three configurations of 0C,IC and TO 
especially in the lower order nodes. It is also observed 
frcn the table that the natural frequencies are clustered 
in the- axi symmetric modes of vibration, especially in the 
lower order nodes, because in the equations (5.16), (5. 17) 
and (5.13) the tern corresponding to C/r is very much 
larger compared to the other two terms that depend on the 
modal number. The loss factors are small and decrease mono- 
tonically with modal number . 

Table 5-2 shows the resonant frequencies and modal 
loss fact ora of the layered shell in the three configura- 
tions of the applied damping treatment with the storage 

♦ Q 2 

modulus of the damping layer changed to = 4.14 x 10 N/m . 
The other parameters are the same as that in Table 5.1. In 
this case also the natural frequencies are clustered near 
the first fundamental frequency. The modal loss factors 
decrease mo not oni cally with modal number . The difference 
between the two results is only in the modal loss factor 
value, which increases by the same order of magnitude as 
the increase in . 

Table 5.3 shows the natural frequencies and loss fac- 
tors for two different values of the damping larger Young's 
modulus and for a thickness of the damping layer ten times 
the thickness of the base layer, for the outside coating 
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of the damping treatment. It is clearly observed from 
the table that with increase in thickness and increase 
in the storage Young's modulus of the damping layer, the 
modal loss factors are considerably increased almost in 
a linear manner both with respect to the damping layer 
thickness, and its Young's modulus. The clustering of 
the natural frequencies around its fundamental frequency 
is seme what reduced with increased values of and 
The modal loss factors again decrease with modal number. 

b) Sandwich Shell : Table 5.4 presents the resonant fre- 
quencies and loss factors of the sandwich, shell for the 
three values of the core shear parameter. In this case 
too, the natural frequencies are closely packed near the 
first resonant frequency. The predominant term in the 
frequency equation (5.31) is the term 

D. 12(1- v 2 ) 



for the first few modal numbers and is independent of 
modal number . This term corresponds to the square of 
the ring frequency of a homogeneous shell having flexu- 
ral rigidity equal to twice the flexural rigidity of a 
shell of face layer material of thickness h^. It is 
because of the independence of this term with respect 
to modal number, that the natural frequencies are clus- 
Ihe modal loss factors turn out to be extremely 


tered. 
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small In the first few modes, essentially due to the same 
reason, the imaginary part of equation (5.31) being negli- 
gibly small compared to the real part. The nodal loss fa- 
ctor increases in all cases of g- with nodal number . As 
has been the case with the infinite and finite sandwich 
plates , it is seen from the table that there exists an 
optimum value g which depends on the modal number for 
maximum modal loss fact or. The natural frequencies in- 
crease monotonically with g , the increase being signi- 
ficant only in the higher modes, for reasons mentioned 
above . 

Table 5-5 gives the resonant frequencies and loss 
factors for different values of the core loss factor p. 

The resonant frequencies do not change significantly with 
respect to p for the particular value of g chosen. The 
modal loss factors increase with p . The increase is very 
nearly linear for modal numbers greater than 5. Thus even 
with increase in p, the modal loss factor in the funda- 
mental mode is not improved signi ficantly due to the pre- 
dominance of the third term in equation (5.31) . 

Table 5.6 shows the variation of resonant frequencies 
and loss factors with respect to the geometric parameter Y . 
The resonant frequencies decrease with increasing Y, the 
Increase becoming less for higher nodal numbers. The Ices 
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factors increase with increase in I, but in the first two 
nodes of the structure the loss factor remains extremely 
small. 

Table 5 .7 shows the cavity resonant frequencies for 
the dimensions of the cavity adopted. 

5.5.2 Noise Reduction and Structural Response 
a) Layered Shell : Figure 5.5 shows the variation of the 
noise reduction inside the cylindrical enclosure for two 
different thickness ratios of the damping layer to the 
face layer, and for a constant value of the damping layer 
Young* s modulus E* = 4.14 x 10 N/m . The excitation con- 
sidered is white noise excitation with spectral density 
function S Q . The noise reduction curves show dips of 
negative values corresponding to most of the cavity re- 
sonant frequencies of Table 5.7 in the case of thickness 
ratio of the damping layer to the elastic face lsyer equal 
to 0.1 . There are many structural resonant frequencies 
and cavity resonant frequencies within the range of 650Hz, 
that it is difficult to describe the coupling action bet- 
ween the structure and the cavity modes. For the case of 
a higher thickness ratio of 10, the noise reduction im- 
proves significantly beyond a frequency of 400 Hz, essen- 
tially due to the increased surface mass density of the 
layered shell. Even though dips occur at the cavity 
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Table 5*7 - Acoustic Resonant Frequencies of a Circular 
Cylindrical Cavity (Axi-symetric Modes ) 


i k f iok 
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3 
1 
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2 

4 

33 

0 

4 

1 

2 

3 

4 


0.0 

109.00 

165.00 

198.29 

201.36 

260.33 

291.99 

330.00 

335.39 
347.84 

382.40 
386.58 
416.48 
440.64 

495 .00 

505.00 
507.07 
534.39 
574.70 
625.50 
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FIG. 5.5 EFFECT OF DAMPING LAYER THICKNESS ON NOISE REDUCTION 
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resonances in this case, the minimum values of noise reduc- 
tion .are in the positive range. Sometimes the increase 
in noise reduction value for the higher thickness ratio 
as compared to the lower thickness ratio is as much as 
40 dB (at 625 Hz). Beyond 625 Hz, there are no dips in 
the noise reduction curves, for the cavity modes beyond 
625 Hz have not been included in the analysis. The very 
low value of the noise reduction for l^/b^ = 0,1 625Hz 

is due to the high degree of coupling between a structural 
mode and cavity mode. 

Figure 5-6 shows the noise reduction curves for the 

same thickness ratios, but for a different value of the 

• Q 2 

Young's modulus of the damping layer E 2 = 4.14 x 10 N/m . 

In this case, the dips at cavity resonances are smaller than 

6 2 

for the corresponding values for 4.14 x 10 N/m . The 

increased thickness of the damping layer does not help to 
improve the noise reduction values for Eg = 4.14 x 10^N/m 2 
unlike the case for B * = 4.14 x 10 6 N/m 2 . The fact that at 
the cavity 1 off resonances' the maximum value of the noise 
reduction values are considerably less for E^=4.14xlO y N/m . 
may be due to the cavity modal impedances having higher 
values because of the acoustoelastic coupling. Because 
of the close proximity of the structural resonant fre- 
quencies and the acoustic cavity frequencies in the frequency 
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range considered, it is very difficult to describe the 
coupling action between the structural motion and the 
cavity pressure with respect to the properties of the 
layer materials. Figures 5.7 and 5.8 respectively show 
the shell response sepectral densities corresponding to 
Figures 5.5 and 5.6. It is observed from Figure 5.7 that 
for /h^ = 0.1 and Ilf = 4.14 x 10^ H/m^, the structural 
response peaks at 475 Hz and at 625 Hz. The peak at 625Hz 
is predominant due to the strong coupling of the struc- 
tural mode and cavity mode which coincide at that frequency. 
For h^/h^ = 10, the response levels are considerably lower 
than for h^/h^ = 0.1 because of high modal loss factors . 
In this case the response peaks at 300 Hz which is con- 
siderably higher than the first natural frequency of the 
structure of 269 Hz. 

From Figure 5.8, it is seen that the structural res- 
ponse has low levels corresponding to hg/h^ = 10.0 as 
compared to hg/b^ = 0.1, which is due to the high modal 
loss factors for the higher thickness ratio. The res- 
ponse peak corresponding to the first structural resonant 
frequencies are predominant especially for h 2 /hj= 0 . 1 . 

The structural response peak in this case is slightly 
shifted to the right due to the coupling action from 
462 Hz to 475 Hz. 


4 14 x 10° Him 1 
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t) Sandwich Shell : Figures 5.9 to 5.11 show the varia- 
tion 0 no ise reduction respectively with respect to the 
core shear parameter , loss factor and geometric parameter, 
for an external white noise excitation. The noise reduc- 
tion has negative values corresponding to the cavity re- 
sonant frequencies for all value of g upto a frequency of 
400 Hz as seen from Figure 5-9. The dip at 275 Hz which 
is neither a cavity frequency nor a structural frequency 
is a resonant frequency of the coupled cavity-structure 
system. Corresponding to this frequency the noise reduc- 
tion value is improved for a higher value of g = 5000. 
Otherwise, there is no significant difference in the noise 
reduction values for the different values of g . Even 
though there are three cavity resonant frequencies in the 
range 400 - 500 Hz, the noise reduction curves do not have 
dips in this range due to the coupling action of the cavity 
and the structural motion. Figure 5.10 shows the variation 
of noise reduction with respect to the core loss factor p. 
Again there is no significant variation in the noise reduc- 
tion curves. The noise reduction curves have dips cor- 
responding to the cavity resonant frequencies. But the 
most predominant dip occurs at 275 Hz which is a resonant 
frequency of the coupled cavity system. 

The insensitive ness of the noise reduction with 
respect to the core loss factor is mainly due to the fact 


IQ, 
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FIG, 5 9 EFFECT OF SHFAR PARAMETER ON NOISE REDUCTION ( SANDWICH SHF 1 1 » 
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that in the lew frequency ranges the modal loss factors 
are very low. 


Pig-are 5.11 shows the variation of the noise reduc- 
tion inside the cylindrical shell with respect to the 
geometric parameter. In this case too, the noise re- 
duction is insensitive to the geometric parameter and the 
differences are only marginal. 

. Pigures 5,12 to 5.14 show the variations in struc- 
tural response with respect to , £ and Y respectively. 

In all the three cases the predominant structural response 
peak occurs at a frequency of 275 Hz corresponding to a 
resonant frequency of the coupled cavity, structure system. 
Surprisingly, the structural response curves are also not 
very sensitive to the variations in the sandwich core para- 
meters . It is due to the low values of the modal loss fac- 
tors obtained in the low frequency ranges because of the 
predominance of the ring frequency term in the natural 
frequency expression of equation (5.31). The constrained 
layer damping treatment in the axisymmetric modes of vi- 
bration of the shell in the lower order structural modes 
is not very effective • 

Figure 5.15 compares the noise reduction of a sandwich 
shell with constrained damping layer treatment to that of 
an equivalent layered shell with unconstrained damping 
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treatment applied as an outside coating. In both the 
cases the surface mass density is kept constant. The 
resonant frequencies aid modal loss factors of the equi- 
valent layered shell is shown in Table 5.8, along with 
those of the corresponding sandwich shell. The modal 
loss factor in the equivalent layered shell is less than 
that for the sandwich shell in all the modes except the 
first mode. From Figure 5.15 it is seen that beyond 
400 Hz the sandwich shell h as slightly improved trans- 
mission characteristics . The dips in the noise reduction 
curves in the case of the sandwich shell are less pronoun- 
ced, beyond that frequency. Figure 5.16 shows the cor- 
responding structural response in the two cases. While 
in the case of the sandwich shell the structural peak 
occurs corresponding to the coupled resonant frequency of 
the cavity-structure system, in the case of the equivalent 
layered shell , the response peaks occur by and large with 
respect to the uncoupled resonant frequencies of the 
structure . 


5.6 CONSLOSIONS 

The problems of the noise transmission of a layered 
shell with unconstrained damping treatment and a sandwich 
shell with constrained damping treatment into a closed 
cylindrical cavity are formulated in an acoustoelastic 





Table 5»8 - Resonant Frequencies and loss Factors 
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5 16 COMPARISON OF STRUCTURAL RESPONSE BETWEEN SANDWICH 
AND EQUIVALENT LAYERED SHELL 
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formulation. Only the axisyrmetric vitiations of the 
shell are considered. The coupled cavity-structure equa- 
tions are solved by a matrix inversion scheme. 

The results show that the applied damping treatment 
is not very effective in the axisynmetric modes of vibra- 
tion of the sandwich shell because of the predominance of 
vibration corresponding to the ring frequency of the com- 
posite shell. Hence the noise reduction and structural 
response are insensitive to the changes in the sandwich 
core parameters. The in-vacuo structural resonant fre- 
quency of the sandwich shell is reduced due to the coupling 
action, as the cavity adds to the inertia of the structure . 

In the case of the layered shell significant improve- 
ments in the noise transmission characteristics of the 
shell are possible by increasing the thickness of the 
damping layer and a suitable choice of the Young's modu- 
lus of the damping layer material. The noise transmis- 
sion characteristics and the structural response do not 
vary with the three configurations of the layer treatment 
00, IC and TO. 
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SUMMARY AID CONCLUSIONS 

This thesis presents mainly a theoretical analysis of 
the noise transmission characteristics of sandwich struc- 
tures with constrained damping layer treatment. The dif- 
ferent problems considered in the analysis are : 

1. Sound transmission through an infinite sandwich 
plate, 

2. Noise transmission through an elastically suppor- 
ted finite sandwich plate into a rectangular cavity 

3. The noise transmission through finite circular cylin- 
drical sandwich shell into a closed cavity . 

The noise transmission through a layered shell with 
unconstrained damping treatment into a cylindrical enclo- 
sure is also studied. The influence of the damping treat- 
ments on the structural response of the composite struc- 
ture is investigated in the analysis. The coupling bet- 
ween the acoustic pressure loading and the structural 
motion is adequately treated in all the examples. 

6.1 INFINITE SANDWICH PLATE 

The problem of the sound transmission of a damped 
sandwich plate of infinite extent is formulated on classical 
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lines. The f emulation takes into account the interaction 
between the structural motion, and the radiated pressures 
on either side of the panel, the incident and reflected 
pressures . 

The damping layer consists of a viscoelastic material 
of relatively low rigidity constrained between two elastic 
face layers . In the analysis the face plates are assumed 
to be of the same material and same thickness . The equa- 
tions of notion of a three layered sandwich plate [6] have 
been reduced to represent the equation of notion of the 
infinite sandwich. 

Expressions for the sound transmission loss of the 
sandwich plate and the coincidence frequency have been 
derived in terms of the sandwich core parameters , the 
face plate properties and the incidence angle of the ex- 
ternal harmonic pressure wave using the dispersion rela- 
tion between flexural wave number and frequency of the 
sandwich plate. 

$ \g - 

Unlike the case of the homogeneous sandwich panel, 
the coincidence frequency for the sandwich panel becomes 
a complex number as a consequence of the flexural wave 
motion in the plate being a spatially damped wave. The 
sound transmission coefficient for the sandwich plate is 
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derived in teres of the specific transmission impedance 
cf the sandwich panel. The specific transmission impe- 
dance of the sandwich panel has both reactive and resis- 
tive components, unlike the case of the homogeneous elas- 
tic panel, whose specific transmission impedance is wholly 
reactive. Because of the damping layer, the transmission 
loss has a positive value at all frequencies. 

A parametric study of the sound transmission loss 
of the sandwich plate is undertaken with respect to the 
core parameters, namely, the shear parameter, core loss 
factor and the geometric parameter, and the incidence angle 
of the incident pressure wave. The following conclusions 
are drawn as a result of the investigation. 

1. The sound transmission loss is very much sensi- 
tive to the variations in the core shear para- 
meter g‘ . The coincidence frequency decreases 
with increasing values of g* . The analysis 
also shows the existence of an optimu m value 
of g* for maximum transmission loss at coin- 
cidence . 

2. The coincidence transmission loss of the sand- 
wich. increases with increase in core loss factor. 

The m axi mu m increase of transmission less with 
core loss feet or is achieved for the optimum 
value cf the shear parameter g* . 
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3 . The transmission loss increases marginally with 
increase in the geometric parameter. 

4. The coincidence frequency for the sandwich panel 
is higher than for an equivalent homogeneous panel 
of the sane surface mass density. 

5. The sound transmission loss decreases with increase 
in incidence angles up to the coincidence frequency. 

The coincidence frequency decreases with increasing 
angles of incidence. 

The study clearly reveals that significant improve- 
ments in the sound insulation characteristics of homoge- 
neous panels can be achieved by constrained layer damping 
treatments mainly due to two reasons . The coincidence 
frequency can be shifted to higher ranges of frequency by 
the sandwich construction. Secondly the coincidence trans- 
mission loss can be increased significantly by a proper 
choice of the core shear parameter and loss factor. 

6.2 FINITE SANDWICH MB 

The noise transmission through a finite sandwich 
plate with constrained damping layer treatment into a 
rectangular cavity is formulated as an acouBt celas ti c 
problem on the lines of Dowell [5]» taking into consi- 
deration the interaction between the internal sound field 
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In the cavity and the structural notion. A direct natrix 
Inversion scheme is proposed and used to solve the coupled 
natrix acoustoelastic equations in the pressure coeffici- 
ents and the structural generalized coordinates. Shis 
method takes advantage of the diagonal nature of some of 
the matrices in a ’partition’ analysis which avoids the 
use of an eigenvalue analysis which nay be a difficult 
proposition in view of the order of the natrix involved 
and the elements in the matrix being complex values . 

The formulation admits elastic boundary conditions 
at the longer edges of the plate representing realistic 
boundary conditions that would obtain in s beet -stringer 
panel construction. The other two edges are assumed to 
be simply supported. The presence of the damping layer 
and the assumption of boundary conditions other than simply 
supported necessitates a ’forced damped normal mode’ ana- 
lysis for the response and noise transmission calculations. 
The damped normal mode analysis is systematically applied 
analogous to the classical normal mode analysis for deter- 
mining the vibroacoustic response. The natural frequen- 
cies, assos$.cated loss factors and damped normal modes of 
the sandwich plate are obtained by an iterative interpola- 
tion method. Expressions for the modal parameters, like 
the generalized mass , generalized force, acoust o-structural 
modal coupling coefficients are derived in terms of the 
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sandwich properties. Only the symmetric structural nodes 
are considered in the analysis in view of symmetric boun- 
dary conditions and uni fern loading on the plate. The 
frequency range of interest in the study has been res- 
tricted to 10C0 Hz. Both harmonic and stationary random 
loading are considered in the analysis. 

An exhaustive parametric study with respect tc the 
sandwich core parameters is undertaken. The following 
conclusions are drawn from the parametric study. 

1. The effect of the boundary conditions on the 
noise reduction is mainly felt only at the fun- 
damental structural resonance** of the panel. 

2. The induced damping improves the noise reduction 
values at structural resonances especially in the 
lower order structural modes. 

3. Negative noise reduction values inside the enclo- 
sure are obtained corresponding to lew structural 
resonant frequencies aid loss factors, and whenever 
the in-vacuo structural resonant frequencies and 
cavity resonant frequencies occur in close proximity. 

4. As in the case of the infinite sandwich plate, 

by a suitable choice of g* , the core shear para- 
meter and the core loss factor p, the modal loss 
fact ora can be increased with consequent improve- 
ment in the noise transmission characteristics cor- 
responding to structural resonant frequencies. 
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5. The structural response reduces significantly with 
increase in the core loss factor p. The structural 
response has maximum values corresponding to its 
first resonant frequency which may be slightly 
altered due to the coupling action of the cavity 
pressure . 

6. The structural response attains higher values also 
corresponding to a coupled natural frequency of the 
cavity-structure system due to close proximity of 
the in-vac uo structural resonant frequency and a 
cavity resonant frequency. 

The results of the vihro-acous t ic response study of 
cavity backed finite sandwich panel are qualitatively si- 
milar to those of finite homogeneous panels , but because 
of the presence of the damping layer the noise transmission 
of the sandwich panel is improved especially at the lower 
order structural modes. 

6.5 FINITE LITERS I) I3D SANDWICH SHELLS 

The problems of noise transmission and structural 
response of a layered circular cylindrical shell with 
unconstrained damping layer treatment and a sandwich 
shell with constrained damping treatment .are considered 
in the analysis in a similar manner as that of a finite 
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sandwich plate. Only the axi -symmetric vibrations of the 
shell and simply supported end conditions are investigated. 

The acoustoelastic problems are formulated in a matrix 
format taking into account the interaction of the cavity 
pressure field and the flexible wall. The cavity pressure 
is expanded in terms of the cavity normal modes (axi- 
symmetric), the radial variations in terms of Bessel func- 
tions of the first kind and the longitudinal variations 
in terms of sine functions. The structural motion is ex- 
panded in terms of sine functions , resulting in terms 
of a coupled matrix differential equations in the pressure 
coefficients and structural generalized coordinates. The 
problem is solved by a matrix inversion procedure similar 
to that adopted for the sandwich plate. It may be mentioned 
here, that, perhaps for the first tine in this work, that 
the noise transmission through layered shells with const- 
rained and unconstrained damping treatments is analysed. 

For the case of the finite shell , extensive para- 
metric studies as in the case of the sandwich plate are 
not carried cut. Moreover, because of the close proxi- 
mity of a number of structural and cavity resonant fre- 
quencies in the frequency range of interest, it becomes 
difficult to explain the coupling action between the ca- 
vity pressure and the structural motion precisely. 
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Damping effectiveness in the case of sandwich shells 
is small xn the lower order axisymmetric modes of vibration. 
In the natural frequency expressions for the sandwich shell, 
a term similar to the ring frequency of homogeneous elastic 
shells, which is independent of the damping layer proper- 
ties and the modal numbers , predominates over the other 
terms. Hence, the natural frequencies are clustered near 
the first resonant frequency and the modal loss factors 
in the first resonant frequency and for subsequent lower 
order modes are very small. 

In the case of the layered shell the modal loss fac- 
tor decreases with modal number, while for the sandwich 
shell the modal loss factor increases with modal number. 

The modal loss factor for the layered shell can be increa- 
sed by increasing the elastic modulus of the damping layer 
and its thickness. However, the modal loss factor asso- 
ciated with the first node of the sandwich shell cannot 
be increased sufficiently as explained earlier. Hence , 
the noise reduction in the case of the sandwich shell is 
quite insensitive to variations in the sandwich core 
parameters . The structural response of the sandwich shell 
is also not affected by the variations in the core para- 
meters. But, the coupled natural frequencies of the struc- 
ture-cavity system are significantly different from the 
in-vacuo structural resonant frequencies for the sandwich 
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shell considered in the analysis. For the layered shell 
however, the shifts in the natural frequencies fr n their 
in— vacuo values are not as much pronounced . 

The noise reduction for the layered shell is sensi- 
tive to both the thickness and the elastic modulus of 
the damping layer, especially so for the non-resonant 
transmission of sound. 

6.4 SCOPS FOR FURTHER RESEARCH 

The problem of noise transmission of sandwich struc- 
tures, presented in this thesis can be extended to include 
irregular geometries of the cavity and arbitrary shapes of 
the sandwich panel. It can be achieved by considering 
finite element formulations both for the sandwich panel 
and the irregular acoustic enclosure . Graggs [68,95*96] 
has used acoustic finite elements and plate finite element 
for the response and noise transmission characteristics 
of flexible homogeneous panels. The finite elements for 
the sandwich panel can be adopted from the works of Lu, 
Killian and Everstine [ 28 , 29 ], Using a combination of 
these finite elements the acoustoelastic problem can be 
solved. 

In the pr obi an formulation it has been consistently 
assumed that the loss factor and the elastic moduli! of 
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the viscoelastic damping layer are frequency and tempe- 
rature independent. But this restriction on the material 
properties of the damping layer can he relaxed in the fi- 
nite element formulation by incorporating the frequency 
variations of the properties of the viscoelastic layer. 

Even though a matrix inversion scheme has been adopted 
for solving the acoustoelastic problem, better understand- 
ing of the precise nature of the coupling between the struc- 
tural motion and the cavity pressure can be had by perfor- 
ming a matrix eigenvalue analysis which will give the coup- 
led natural frequencies of the cavity-structure system and 
the corresponding loss factors. The method given by 
Meirovich [93] can be extended to handle matrices having 
complex elements . 

In the analysis of the sound transmission of finite 
plates and shells the frequency range of interest has been 
restricted to 1000 Hz. The frequency range of interest 
can.be extended to higher frequencies. In that case the 
equations of motion of the sandwich plate and shell adop- 
ted have to be modified to include the effects of shear 
deflection and rotary inertia. For very high frequency 
ranges the a ccust ©elasticity analysis presented in the 
thesis will not be suitable as one may have to contend 
with a large number of structural and cavity modes. The 
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statistical energy analysis (SEA) may be used for the pre- 
diction of response and noise transmission statistics in 
such cases. The SEA parameters, like the modal density 
and the radiation impedance of the sandwich panels can be 
derived analytically from the appropriate equations of motion 
of the sandwich structures. However, sufficient care has 
to be exercised in the energy flow calculations, as the 
usual assumption of weak coupling .between the resonant 
modes in the frequency bands is not valid, due to the 
large modal damping induced by the constrained damping 
layer with significant intermodal coupling. Recently, 

Chand iramani [97] and Smith [98] have presented SEA models 
where the transition from weak to strong coupling between 
the modes are analysed. Making use of these models it may 
be possible to solve the couple d acoustoelastic problem of 
damped sandwich structures in the high frequency ranges. 

In the case of the layered and sandwich shells, only 
the axisymmetric vibrations are considered. The formula- 
tion can be extended to consider the circumferential lobar 
modes of the shell. 
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APPSIJDH A 


BOUNDARY CONDITIONS AND JUMPED NORMAL MCB3S 

The symmetric modes are obtained by applying the 
boundary conditions at one edge only, namely at x = a/2. 

(a) The boundary conditions for edges with elastic rotational 

restraint and no transverse displacement and unriveted 
edges are , 

3 A r 

2 A_ cos — - = 0 (A.l) 

r=l L 2 

Ia + ^ )[a ' + + g(i+y)/(A ^ + g)]coD 4 s 

+ p +2 SV ) A r s i» =0 (A.2) 

2 Aj»[ ( Ay + P n )/( + P Q +g5] cos - ss 0 (A. 3) 

r=l 2 

The above three boundary conditions are due to zero dis- 
placement, moment equilibrium between stiffener and plate 
and no axial stress in either of the fhce plates respec- 
tively. Ep s E s p /(D^a •*), is the torsion bending stiff- 
ness parameter and K sv = /(D^. a) is the St .Tenant torsion 

stiffness parameter of the stringer cross-section, B g .and 

G- are the elastic and shear moduli! of the stringer 

a 

material, P , the torsion bending constant and J is 
St .Tenant torsion constant of stiffener cross-section. 
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(t>) The first two boundary conditions for edges with 
elastic rotational restraint and no transverse displace- 
ment, but with riveted edges are the same as given by 
equations (A.l) and (A.2). The last boundary condition 
arises because of the prevention of shear strain in the 
core due to riveting. This is given by 


5 y A r ( A 2 + P 2 }] A r 

2 — 2__ sin — ^ = 0 

1=1 [ A 2 + ? 2 + g ] 2 


(A. 4} 


{ c) The boundary condition for transverse elastic con- 
straint and no rotation are 


^ ij.ApMi.-J: - o 


3 y A r ( * 2 + 

2 <5 5 

r*l + P n + s 1 



(A.5) 


(A. 4) 


3 

* *r 
r*l 


[ + (2- tfBJ] . 


[ 


4 + 4 + 8 ( h!l ]3in 


Ay + 


p2 


Kf cos -jr - ] = 


+ g 


( A*6) 

where = (BI) S / (a D^) is the flexural stiffness 
parameter. Equation (A.6) represents the equilibrium 
of shear forces between the stiffener and the plate. 


0 
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(BI) S is the flexural rigidity of the stiffener. 

(d) The boundary conditions for the fixed-fixed edges 

at x = + a/2, are 


3 A 

£ A cos — — 
r=l r 2 


= 0 


(A.l) 


5 * *r 

2 A jk sin — 

r=l * r 2 


0 


3 V *r< + 4)] 

r=l [ Ar + E n + S ^ 


sin 


= 0 


(A.5) 


U. 4) 


In each of the above cases the farced damped normal modes 
can be expressed as 




(coS A lmn’ + M ros A 2mn ?+ ScoS A 3mn f }si; 


(A.7) 

where, A lTnT1 , A 2mn and are the complex roots of 

the bicubic equation (2.15) corresponding to the mn mode 
(The other three roots are just the negative of A 


Inn* 2nn 

and A^ mn ). The quantities K and 11 corresponding to the 
different edge conditions are obtained as. 


(a) So transverse displacement and elastic rotational con- 


straint with no rivets : 

A 

( R^-R^ ) cos 


v lmm 

T 


(R^-Rg) 


cos 


frignn ' 

“2 


, H 


<V®2 ) 


cos 


_lmn 


(J^-R^jcos 


' Mfaa r 


(A.9) 


2 
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(b) Ho transverse displacement 

constraint with rivets, 

A, 


rivets, 
sln “ A 3 mi ®3 008 


and. elastic rotational 


tan-fSS 

3nn V 0S A 

SmA cos -f®* 8 ” -V 2 - 


; an h2ss _ A o sln h|SS 

2 2 fcinn ^ 


•*Wi sin -f® 


(A. 10) 


'2-A^ ^ tan A|3n . ^ ^ sin A|sn 


(c) Transverse elastic constraint and no rotation 


(A.l 1 ) 


M = 


H = 


Sin -J|5 R 3 A lmn^ 

sin (^3 ^2mn ~ ? *2 *3mn^ 

sin K lmn - ^ A 2an ) 

sin ^2mn “ ^3mn^ 


(A. 12 ) 


£ Jt£* #i*3 ^ 


(d) Fixed edge conditions 

X lnn 8in (S 1 - R 3 ) 


M— 


A 2mn 6in -V ES(R 3- R 2 > 


sm -^aa - Rg) 


‘Iran 


Hs* 


4m 


(®2 Ej ) 


(A. 14 } 


(A. 15 ) 


^3mn s ^ n 


2 
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In the above equations 

2 2 
A~_ + p* 

V s ££2 S f r = 1,2,3 (A.16) 

A + P + ^ 

rmn n 6 

and *t-4 x p + W is the equivalent torsional 

stiffness parameter. 

It should be noted that the results for the simply 
supported edge conditions can be obtained by letting 
M = N = 0 and A lmy( = cm, {m = 1,3,5..) in equation (A.7). 



APPENDIX B 


GENERALIZED MASS , GENERALIZED FORCE, ACCEP- 
TANCE FUNCTIONS AND COUPLING COEFFICIENTS 


The egressions for the generalized mass and the 
generalized force from equations (4.15) (2,31) and (2.32), 
corresponding to ihe forced damped normal modes are. 


% A 

pah ( 1+M^+N 2 ) sin( — — - 4 ^SS ) 


*W 


f- 


2 i 


+2M [■ 


A lmn + A 2m n 


sln( A laa 


A lmn “ A 2mn 


i sim 
] +2N[ 


A lmn + A 3mn \ 

* ) 


A +■ A, 

Imn 3nm 


sm( ±wg lAa aU) 8ln( 
]+2MN [ — 2mn - - 


A 2mn + A 3mn 

— — * ' g -- *— 


’) 


v lmn 


K 3mn 


Bia( A2nn 2 ~ x ?° £ ) 
A 2mn” A 3mn 


-14- 


sin A 


2 A 


Imn + jj2 


3mn 
sin A 


Imn 


TT 

*•» 


2 an 


2mn 


-sin A 


3mn 


2 *3mn * 
Pmn (t) = p o elp [ wt hlm 


(B.1) 


(B.2) 
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where , 


2ab 5 a y sin — ~^c 

*■» - ^ £ - X — j ™ 

1=1 < 4n - 


I - Foos-^pstaf) 


(B.5) 


where, = 1, c 2 = M and Oj = N and F = —-sin© 


Note for n even, y wt , = 0 and along the x direction only 
the symmetric modes are considered. 

Por r = 1, from equation (2.36), the expression for (to ) 

is independent of frequency and is given by 


X = Y Y* 

mnrs 'mn r rs 


(B.4) 


The acoustic structural coupling coefficient given in 
equation (4.21) is obtained from the damped normal modes 


^mnijT" 


2abn (l-cosnrt cos jit) ^ X- 

£ a ■ 

% (n -j ) r=l r 


^ *K 

' ran 


2 , i* i 

cos (vp - ) 


(B.5) 


Because only the synmetric structural modes are admissible, 
the contributions to the acoustic-structural modal coupling 
coefficients are only from the even order acoustic modes 
in the x and y directions. 



